Conf Univ Dr Dana Constantinescu Ecua􀄠ii Diferen􀄠iale Elemente teoretice şi aplica􀄠ii Editura Universitaria, CUPRINS Prefa􀄠ă Considera􀄠ii generale Introducere Noţiuni fundamentale Exerciţii propuse Ecua􀄠ii diferen􀄠iale de ordinul I Ecuaţii cu variabile separabile Ecuaţii liniare Ecuaţii cu diferenţială totală Ecuaţii reductibile la ecuaţii fundamentale Exerciţii propuse Aplica􀄠ii ale ecua􀄠iilor diferen􀄠iale de ordinul I Suprafaţa de echilibru a unui lichid în rotaţie (integrare directă) Evolutia unei populaţii (ecuaţie cu variabile separabile) Variaţia presiunii atmosferice în raport cu altitudinea (ecuaţie cu variabile separabile) Căderea liberă (ecuaţie cu variabile separabile) Descărcarea unui condensator într-o resistenţă (ecuaţie liniară) Incãrcarea unui condensator printr-o rezistenţã în prezenţa unei surse de curent continuu (ecuaţie liniară) Transformarea energiei electrice in căldură (ecuaţie liniară) Formula fundamentală a curentului alternativ (ecuaţie liniară) Oglinda parabolică (ecuaţie omogenă) Exerciţii propuse Ecua􀄠ii diferen􀄠iale de ordin superior Ecuaţii liniare Ecuaţii liniare cu coeficienţi constanţi Ecuaţii liniare cu coeficienţi variabili Ecuaţii incomplete Exerciţii propuse Sisteme de ecua􀄠ii diferen􀄠iale de ordinul I Consideraţii generale Sisteme liniare omogene cu coeficienţi constanţi Sisteme liniare neomogene cu coeficienţi constanţi Exerciţii propuse Elemente de calcul opera􀄠ional şi aplica􀄠ii în teoria ecua􀄠iilor diferen􀄠iale Transformata Laplace Transformata Laplace inversă Calcul operaţional Rezolvarea ecuaţiilor liniare cu coeficienţi constanţi Rezolvarea sistemelor de ecuaţiil liniare cu coeficienţi constanţi Exerciţii propuse Aplica􀄠ii ale ecua􀄠iilor diferen􀄠iale de ordin superior Oscilaţii armonice Mişcarea unui pendul (ecuaţie liniarǎ) Calculul perioadei unui circuit oscilant (ecuaţie liniarǎ) Oscilaţii ale unei coloane de lichid (ecuaţie liniarǎ) Propagarea cǎldurii într-o barǎ (ecuaţie liniarǎ) Ecuaţii de mişcare (ecuaţii liniare) Determinarea coeficientului de frecare (ec neliniară) Determinarea ecuaţiei unei curbe (ecuaţie neliniară) BIBLIOGRAFIE ………………………………………………… PREFA􀄟􀓼 Materialul de faţă se doreşte a fi o introducere în studiul ecuaţiilor diferenţiale şi al aplicaţiilor lor Lucrarea se adresează în primul rând studenţilor facultăţilor tehnice, precum si tuturor celor ce doresc să folosească noţiuni fundamentale privind teoria ecuaţiilor diferentiale în aplicaţii practice Scopul lucrării este prezentarea clară şi precisă a noţiunilor şi rezultatelor de bază privind rezolvarea analitică a unor tipuri importante de ecuaţii diferentiale ordinare, descrierea şi exemplificarea principalelor metode de rezolvare a problemelor, precum şi folosirea acestora pentru modelarea matematică a unor probleme practice Lucrarea are un pronunţat caracter metodic Materialul prezentat este organizat in capitole care acoperă cunoştinţele aferente temei „Ecuaţii diferenţiale ordinare”, temă predată şi seminarizată la disciplina fundamentală Matematici Speciale, prevăzută la unele facultăţi tehnice în noul plan de învăţământ, în semestrul II, anul I Din această perspectivă lucrarea este un eficient suport pentru pregătirea seminariilor, temelor de casă şi examenelor Consideraţii generale asupra ecuaţiilor diferenţiale ordinare sunt prezentate în Capitolul Principalele tipuri de ecuaţii diferenţiale de ordinul I sunt analizate în Capitolul , iar rezolvarea ecuaţiilor de ordin superior reprezintă subiectul Capitolului Legătura dintre acestea şi sistemele de ecuaţii de ordinul I este analizată în Capitolul , unde sunt prezentate şi metodele de rezolvare a sistemelor de ecuaţii liniare de ordinul I cu coeficienţi constanţi Capitolul arată în ce măsură calculul operaţional poate fi folosit pentru rezolvarea (analitică) mai simplă a unor categorii speciale de ecuaţii diferenţiale, cum sunt ecuaţiile şi sistemele de ecuaţii liniare cu coeficienţi constanţi, ce stau la baza studiului sistemelor dinamice O atenţie deosebită este acordată interpretării noţiunilor introduse şi prezentării unor situaţii în care ele sunt folosite Capitolul şi Capitolul sunt dedicate exclusiv prezentării unor aplicaţii importante ale ecuaţiilor diferenţiale ordinare, care reprezintă un element extrem de important în înţelegerea unor fenomene foarte variate (ce apar în fizică, mecanică, inginerie electrică, ecologie, etc ) Descrierea suprafeţei de echilibru a unui lichid în rotaţie, variaţia presiunii atmosferice în raport cu altitudinea, studiul ecuaţiilor de mişcare a corpurilor, a unor fenomene legate de circuitele electrice (descărcarea unui condensator într-o rezistenţă, încărcarea unui condensator printr-o rezistenţă în prezenţa unei surse de curent continuu, calculul perioadei unui circuit oscilant), descrierea propagării căldurii într-o bară, sunt doar cteva dintre aplicaţiile analizate Aceste aplicaţii au fost alese din perspectiva pregătirii studenţilor cărora materialul li se adresează cu precădere Materialul este conceput într-o manieră accesibilă şi sistematică, în aşa fel încât să poată însoţi primii paşi în descoperirea fascinantei lumi a ecuaţiilor diferenţiale Demosntraţiile teoremelor sunt reduse la strictul necesar unei prezentări riguroase, fără încarce discursul matematic Sperăm ca folosirea acestui material să fie doar începutul studiului în domeniul inepuizabil al ecuaţiilor diferenţiale şi să trezească interesul pentru matematicile aplicate, îmbinare de rigurozitate, ingeniozitate şi pragmatism Mulţumim tuturor celor care, prin sugestii şi completări, au contribuit la realizarea acestei lucrări, în special referenţilor pentru observaţiile şi aprecierile făcute cu ocazia analizei manuscrisului CAPITOLUL Considera􀄠ii generale Introducere Studiul ecuatiilor diferenţiale formeazã obiectul unui capitol foarte important al matematicii, atât datoritã rezultatelor teoretice deosebit de interesante cât şi pentru cã ele au nenumãrate aplicaţii în cele mai diverse domenii Ceea ce deosebeşte o ecuaţie diferenţialã de o ecuaţie algebricã este faptul cã necunoscuta nu este un numãr ci o funcţie care satisface o anumiã egalitate şi care trebuie determinatǎ Multe fenomene sunt descrise cu ajutorul ecuaţiilor diferenţiale obţinute prin metoda cunoscutã sub numele de “metoda diferenţialelor” Aceasta constã în înlocuirea unor relaţii ce apar între creşterile infinit de mici ale unor cantitãţi care variazã în timp prin relaţii între diferenţialele (derivatele) lor Spre exemplu viteza instantanee ( ) v t de deplasare a unui mobil care la momentul t a parcurs distanţa s (t ) este ( ) ( ) ( ) ( ) lim ' t t s t s t v t s t → t t − = = − La rândul sãu acceleraţia corpului la momentul t este ( ) ( ) ( ) ( ) ( ) lim ' '' t t v t v t a t v t s t → t t − = = = − In relaţiile ce descriu mişcarea viteza se va considera v(t ) = s '(t ) şi a(t ) = v'(t ) = s''(t ) Exemplu : Mişcarea unui corp sub ac􀄠iunea greutã􀄠ii sale şi întâmpinând o rezisten􀄠ã a aerului propor􀄠ionalã cu viteza sa (acest caz corespunde vitezelor mici) poate fi descrisã cu ajutorul unei ecua􀄠ii diferen􀄠iale Se noteazã v(t ) viteza instantanee a corpului la momentul de timp t > Rezistenţa aerului va fi R(t ) = kv(t ) Legea fundamentalã a mecanicii ( F = ma ur r ) conduce la relaţia mg − kv(t ) = mv '(t ) care reprezintã o ecuaţie diferenţialã cu necunoscuta v = v(t ) Pentru a determina viteza instantanee a corpului trebuie rezolvatã aceastã ecuaţie Probleme fundamentale în teoria ecuaţiilor (în general) sunt determinarea soluţiilor lor sau aproximarea acestor soluţii dacã determinarea analiticã nu este posibilã Teoria ecuaţiilor diferenţiale are mai multe ramuri: - teoria cantitativã se ocupã de rezolvarea analiticã a ecuaţiilor Sunt precizate tipurile de ecuaţii ale cǎror soluţii se pot obţine analitic şi tehnicile de rezolvare a lor - teoria calitativã încearcã sã deducã proprietãţile soluţiilor, chiar dacã expresia lor analiticã nu poate fi cunoscutã - aplicarea metodelelor numerice pentru aproximarea solu􀄠iilor Scopul acestui capitol este prezentarea celor mai importante elemente ale teoriei cantitative a ecuaţiilor diferenţiale No􀄠iuni fundamentale Defini􀄠ia Se numeşte ecua􀄠ie diferen􀄠ialã cu variabila independentǎ x , şi funcţia necunoscutã y = y (x) o egalitate de forma ( , ( ), '( ), , ( ) ( )) n F x y x y x y x = ( ) unde F :D ⊂ Rn+ →R este o functie datã, continuã pe domeniul sãu de definiţie, iar ', '', ( ) n y y y sunt derivatele lui y Dacǎ ecuaţia ( ) este scrisã sub forma y (n) = f (x, y(x), y' (x), , y (n− ) (x)) ( ’) se spune cã are formã explicitã Ecuaţia diferenţialã are ordinul “n” dacã derivata de ordin maxim care apare in ecuaţie este (n) y Exemple: ) x + y (x) = nu este ecuaţie diferenţialã, pentru cã derivatele funcţiei necunoscute “y” nu apar in ecuatie, dar x + ( y' ' (x)) = este o ecuaţie diferenţialã de ordinul cu funcţia necunoscutã “y” şi variabila independentã “x” )mv' (t ) + kv(t ) − mg = este o ecuaţie diferenţialã de ordinul I cu necunoscuta “v” şi variabila independentã “t” Ea descrie mişcarea unui corp sub acţiunea greutãţii sale şi întâmpinând o rezistenţã a aerului proporţionalã cu viteza sa (funcţia necunoscutã, v , este viteza corpului) ) ( ) ( ) C R q t q t ⋅ ' = − este o ecuaţie diferenţialã de ordinul I cu necunoscuta “q” şi variabila independentã “t” Ea descrie procesul de descãrcare al unui condensator de capacitate C într-o rezistenta R (funcţia necunoscutã “q” reprezintã sarcina electricã) ) x ⋅ ln x − y'+ y"− xy + x y' ' '= reprezintã o ecuaţie diferenţialã de ordinul cu necunoscuta “y” şi variabila independentã ‘x” ) (x + y + )dx + (x − y + )dy = reprezintã o ecuaţie diferenţialã de ordinul I pentru cã apar notatiile « dx » şi « dy » asociate formal cu derivatele de ordinul I Necunoscuta problemei trebuie precizatã : dacã folosim notaţia y'= dv / dx atunci ecuaţia se scrie sub forma (x + y + ) + (x − y + )y'= şi necunoscuta ecuaţiei este y, dar dacã vom considera cã x'= dx / dy ecuaţia se scrie sub forma (x + y + )x'+(x − y + )= şi necunoscuta ecuaţiei este « x » Defini􀄠ia Se numeşte solu􀄠ie a ecua􀄠iei diferen􀄠iale ( ) pe intervalul I ⊂ R orice funcţie φ : I →R , derivabilã de n ori pe I , care verificã ecuaţia, adicã pentru orice x∈I are loc egalitatea ( , ( ), '( ), , ( ) ( )) n F x φ x φ x φ x = Existã trei tipuri de solutii: • Solu􀄠ia generalã a ecuaţiei ( ) este soluţia care depinde de x şi de n constante arbitrare n C ,C , ,C (exact atâtea cât este ordinul ecuaţiei), adicã este de forma ( ) , , , , n y =φ x C C C Aceasta este forma explicitã a soluţiei pentru cã se precizeazã modul în care funcţia necunoscutã y depinde de variabila independentã x Uneori soluţia generalã este prezentatã în formã implicitã (integrala generalã a ecuaţiei), adicǎ ( ) , , , , n Ω x y C C = Soluţia generalã se poate obţine şi sub formã parametricã : ( ) ( ) , , , , , , , n n x = f t C C y = g t C C • Orice soluţie care se obţine din soluţia generalã pentru anumite valori particulare ale constantelor se numeşte solu􀄠ie particularã • Soluţiile ecuaţiei care nu se pot obţine prin acest procedeu din soluţia generalã se numesc solu􀄠ii singulare In probleme practice, alãturi de ecuaţia diferenţialã trebuie considerate şi condiţii iniţiale ( ) ( ) ( ) ( ) , ' , , n n y x y y x y y x y − = = = − ( ) Ecuaţia ( ) împreunǎ cu condiţiile iniţiale ( ) formeazǎ o problemã Cauchy Soluţia unei probleme Cauchy ( )+( ) se obţine impunând condiţiile iniţiale ( ) soluţiei generale a ecuaţiei ( ) Exemple y ' = x este o ecuaţie diferenţialã de ordinul I cu necunoscuta y Soluţia sa generalã este y : R→R, y (x) = x + C deoarece verifica ecuaţia Ea depinde de o singurã constantã y (x) = x + , y (x) = x − sunt soluţii particulare ale ecuaţiei pentru cã au fost obţinute din soluţia generalã pentru C = , respectiv C = − Existã o infinitate de soluţii particulare ale ecuaţiei y ' = − y este o ecuaţie diferenţialã de ordinul I cu necunoscuta y si variabila independentã x Funcţia : , , ( ) sin ( ) y C C R y x x C ⎡ π π ⎤ ⎢ − + ⎥→ = − ⎣ ⎦ reprezintã soluţia generalã a ecuaţiei Funcţia : , , ( ) sin y Ry x x ⎡ π π ⎤ ⎢− ⎥→ = ⎣ ⎦ este soluţie particularã (obţinutã din soluţia generalã pentru C = ) Funcţia : [ , ] , ( ) sin cos y R y x x x π π → = ⎛⎜ − ⎞⎟ = − ⎝ ⎠ este soluţie particularã (obţinutã din soluţia generalã pentru C π = ) Alte soluţii particulare se pot obţine în acelaşi mod, pentru fiecare domeniul de definiţie fiind altul Ecuaţia admite soluţiile singulare ( ) y : R→R, y x = şi ( ) y : R→R, y x = − ( ) ( ) ( ) ( ) y − y + y − y = este o ecuaţie diferenţialã de ordinul Soluţia sa generalã este ( ) y : R→R, y x = C + C x + C e x + C cos x + C sin x y (x) = x, y (x) = − sin x sunt exemple de soluţii particulare Problema Cauchy ( ) ( ) ( ) ( ) ''' '' ' , ' , '' , ''' yIV y y y y y y y ⎧ − + − + = ⎪⎨ ⎩⎪ = = = = are soluţia y (x) = xe x + cos x Aceastã soluţie se obţine din soluţia generalã a ecuaţiei diferenţiale, anume ( ) y x = C e x + C xe x + C cos x + C sin x determinând constantele din sistemul ( ) ( ) ( ) ( ) ' '' ''' y C C y C C C y C C C y C C C ⎧ = + = ⎪ = + + = ⎪⎨ = + − = ⎪⎪ ⎩ = + − = Soluţia sistemului este , , , C = C = C = C = deci soluţia problemei Cauchy este y(x) = xe x + cos x Soluţia generalǎ a ecuaţiei y ⋅ ln y + (x − ln y)⋅ y'= , satisface relaţia ⋅ x ⋅ ln y = ln y + C Aceasta este forma implicitǎ a soluţiei In adevǎr, prin derivarea relaţiei anterioare obţinem ' ' ln ln y y y y y x ⋅ y + ⋅ = ⋅ ⋅ ⋅ , de unde rezultǎ y ⋅ ln y + (x − ln y)⋅ y'= adicǎ faptul cǎ funcţia y satisface ecuaţia diferenţialǎ şi deci este soluţia sa generalǎ (deoarece depinde de o constantǎ) Determinarea formei sale explicite este mai dificilǎ Soluţia generalǎ a ecuaţiei ( x + y + y )+ ( xy + x)y'= satisface relaţia x + x y + x y = C Derivând relaţia anterioarǎ obţinem x + x y + x yy'+ x y + x y'= adicǎ x ⋅[( x + y + y )+ ( xy + x)y']= ceea ce aratǎ cǎ funcţia y satisface ecuaţia (variabila independentǎ x trebuie sǎ şi ia valori nenule deci primul factor al produsului poate fi considerat nenul) O problemã importantã in teoria ecuaţiilor diferenţiale este determinarea soluţiei generale a unei ecuaţii diferenţiale date Acest lucru este posibil numai pentru un numãr restrâns de ecuaţii Unele din aceste cazuri sunt prezentate în paragrafele ce urmeazã Exerci􀄠ii propuse S􀓽 se precizeze dac􀓽 func􀄠ia y = y(x) este solu􀄠ie a ecua􀄠iei diferen􀄠iale (sau a problemei Cauchy) în urm􀓽toarele cazuri Pentru fiecare solu􀄠ie s􀓽 se precizeze domeniul s􀓽u de defini􀄠ie a) y + xy − x y'= y(x) = −x /(C + ln x) b) y'+ xy = xe−x ( ) ( ) x y x = x +C e− c) y'−y = xy y(x)= /( − x +C ⋅ e−x ) d) ( ) ( ) ⎩ ⎨ ⎧ = = − + = , ' '' ' y y y y y y(x) = ex + e x e) y''− y' y = ⋅ ex y(x) C e x C e x x = ⋅ + ⋅ + f) ( ) ( ) ( ) ( ) ⎪⎩ ⎪⎨ ⎧ = = + − + + = , ' '' ' y y x y x y y y(x) = ( x + )⋅ ( x + ) S􀓽 se arate c􀓽 solu􀄠ia general􀓽 a ecua􀄠iei (x − y + )y'+x + y + = , scris􀓽 sub form􀓽 implicit􀓽 este C y x y x y x + + ⋅ + ⋅ − = S􀓽 se arate c􀓽 solu􀄠ia problemei Cauchy x + y ⋅ y'= , y( ) = satisface rela􀄠ia x + y = S􀓽 se precizeze forma parametric􀓽 a solu􀄠iei Solu􀄠ia general􀓽 a ecua􀄠iei y'= − y este y":[C −π / ,C +π / ]→R , y(x) = sin(x −C) S􀓽 se scrie solu􀄠ia problemei Cauchy y'= − y , y( π / ) = / Rezolvare: Din sin( π / −C) = / rezultǎ C =π / deci soluţia este y :[ ,π ]→R, y(x) = sin(x −π / ) = −cos x CAPITOLUL Ecua􀄠ii diferen􀄠iale de ordinul I Ecuaţiile diferenţiale de ordinul I au forma F (x, y, y ') = Cel mai adesea ele sunt scrise în formã explicitã y'= f (x, y) Soluţia lor generalã depinde de o singurã constantã Nu orice ecuaţie diferenţialã de ordinul I poate fi rezolvatã analitic Din punctul de vedere al rezolvãrii analitice existã douã categorii importante de ecuaţii : - ecuaţii fundamentale (ecuaţiile cu variabile separabile, ecuaţiile liniare, ecuaţii cu diferentiale totale) - ecuaţii reductibile la ecuaţii fundamentale (ecuaţii omogene si reductibile la ecuaţii omogene, ecuaţii care admit factor integrant, ecuaţii de tip Bernoulli, de tip Riccati, de tip Lagrange, de tip Clairaut etc) Este foarte importantã cunoaşterea algoritmului de rezolvare a ecuaţiilor fundamentale precum şi a metodelor de reducere a celorlalte ecuaţii la ecuaţii fundamentale Ecua􀄠ii cu variabile separabile Forma generalã a ecuaţiei este y ' = f (x)⋅ g ( y) ( ) unde f , g sunt funcţii reale date, continue pe domeniul lor de definiţie Soluţiile ecuaţiei g ( y) = sunt soluţii, de obicei singulare, ale ecuaţiei Dacã g ( y) ≠ rezolvarea constã in separarea variabilelor urmatã de integrare Metoda de rezolvare: - se rezolvã ecuaţia g ( y) = cu soluţiile , , , k y y y - se scriu soluţiile singulare ale ecuaţiei ( ) ( ) ( ) , , , k y x = y y x = y y x = y Domeniul lor de definiţie este domeniul de definiţie al funcţiei f - se scrie ecuaţia sub forma ( ) f (x) g y y = ' (ceea ce este posibil pentru g(y) ≠ ) şi se obţine integrala generalã a ecuaţiei : ( ) ( ) dy f x dx C g y ∫ = ∫ + , adicã forma implicitã a soluţiei - din integrala generalã se calculeazã (dacã este posibil) y şi se obţine forma explicitã a soluţiei Observa􀄠ie : Soluţia particularã a ecuaţiei ( ) care îndeplineşte condiţia iniţialã ( ) y x = y este datã de ( ) ( ) y x y x ds f t dt g s ∫ = ∫ Ea se poate obţine din soluţia explicitã impunând condiţia iniţialǎ O formã particularã a ecuaţiei cu variabile separabile este y ' = f (x) Soluţia generalã a acestei ecuaţii este y (x) = ∫ f (x)dx Exemple : Sã se rezolve y ' = x + sin x Soluţia generalã este ( ) ( ) sin cos x y x = ∫ x + x dx = − x + C ' x y x = + Soluţia generalã este ( ) ( ) ln x y x dx x C x = = + + + ∫ ' y y x = − In acest caz ( ) f x x = − şi g ( y) = y , deci ecuaţia g ( y) = are soluţia y = şi funcţia y : R −{ }→R, y (x) = s s este soluţie singularã a ecuaţiei Dacã y ≠ ecuaţia devine y ' y x = − şi integrala ei generalã este dy dx y x ∫ = ∫− Rezultã ln | | ln ln ln ln notatie y C C C x x x = + = + = , unde C > este o constantǎ arbitrarǎ Rezultǎ x C y = ± , C > In acest caz soluţia generalǎ se scrie sub forma : { } , ( ) C y R R y x x − → = , C ≠ Inlocuind C = în soluţia generalǎ se obţine soluţia singularǎ s y Aceasta nu este însǎ o soluţie particularǎ deoarece valoarea C = nu este acceptabilǎ în cadrul soluţiei generale Ecua􀄠ii liniare Forma generalã a ecuaţiei liniare este y ' = P(x)⋅ y +Q(x) ( ) unde P,Q: I →R sunt funcţii date, continue pe domeniul de definiţie Aceastã ecuaţie se rezolvã prin metoda variaţiei constantei Metoda de rezolvare (metoda variaţiei constantei) - se rezolvã ecuaţia omogenã y ' = P(x)⋅ y care este o ecuaţie cu variabile separabile şi se obţine soluţia nenulã ( ) y C e C f (x) P x dx notatie x = ⋅ a = ⋅ ∫ - se considerã constanta C ca fiind funcţie de x , adicã se scrie y (x) = C(x)⋅ f (x) - se calculeazã y '(x) = C'(x) f (x) + C(x) f '(x) şi se introduce in ecuaţia ( ) ; termenii care conţin pe C(x) se reduc şi se obţine o ecuaţie mai simplã de forma C'(x) = g (x) - se rezolvã ecuaţia C'(x) = g (x) şi se obţine soluţia C(x) = ∫ g (x)dx + K - se introduce expresia lui C(x) în y (x) = C(x) f (x) şi se obţine forma explicitã a soluţiei ecuaţiei ( ) Observa􀄠ie : Forma explicitã a soluţiei ecuaţiei ( ) , pentru x fixat, este ( ) ( ) ( ) ( ) x x x s P t dt P t dt x x y x K Q s e ds e ⎡ ⎛ ⎞ ⎤ = ⎢ + ⎜ ⎟ ⎥ ⋅ ⎢ ⎜ ⎟ ⎥ ⎣ ⎝ ⎠ ⎦ ∫ ∫ ∫ Aceastã expresie se obţine folosind algoritmul anterior dar e dificil de memorat şi de aceea se recomandã folosirea algoritmului pentru rezolvarea fiecãrei ecuaţii Exemplu : Sã se rezolve problema Cauchy ( ) ' sin / y y ctgx x x y π a = ⋅ + ⋅ ⎧⎪⎨ ⎩⎪ = Funcţia ctgx nu este definitã in punctele nπ , n∈N Din cauza condiţiei iniţiale se va cãuta soluţia generalã a ecuaţiei pe intervalul ( π ) Ecuaţia omogenã y ' = y ⋅ ctgx are integrala generalã cos sin dy x dx y x ∫ = ∫ Rezultã ( ) ln | y |= ln | sin x | +C = ln C | sin x | care dã soluţia y (x) = Csin x Se aplicã variaţia constantei, adicã se considerã y (x) = C(x)sin x Introducând y '(x) = C'(x)sin x + C(x)cos x în ecuaţia neomogenã obţinem ( ) ( ) ( ) cos ' sin cos sin sin sin x C x x C x x C x x x x x + = + Termenii conţinând factorul C(x) se reduc şi se obţine ecuaţia C'(x) = x cu soluţia C(x) = x + K Introducând aceastã expresie în forma lui y (x) obţinem soluţia generalã a ecuaţiei, anume y : ( ,π )→R, y (x) = (x + K)sin x unde K ∈R este o constantã arbitrarã Din condiţia y (π / ) = a rezultã sin K a π π ⎛ ⎞+ ⎛⎞ = ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ ⎝ ⎠ , adicã K a π = − Deci soluţia problemei Cauchy este ( ) ( ) : , , sin y Ry x x a x π π ⎛ ⎞ → = ⎜ + − ⎟ ⎝ ⎠ Ecua􀄠ii cu diferen􀄠ial􀓽 total􀓽 Forma generalã a ecuaţiei este P(x, y)dx +Q(x, y)dy = ( ) unde P,Q sunt funcţii date, de clasã C pe domeniul D ⊂ R şi satisfac relaţia ( ) (x y) x Q x y y P , , ∂ = ∂ ∂ ∂ pentru orice (x, y)∈D Rezolvarea ecuaţiei se bazeazã pe faptul cã existã funcţii de forma ( ) ( ) ( ) , , , x y x y U x y = ∫ P t y dt + ∫ Q x t dt astfel încât ( ) ( ) ( ) , , , x y dU = P x y dx +Q x y dy Spunem în acest caz cã ecuaţia are diferenţialã totalã Ea se scrie sub forma ( ) , = x y dU , deci soluţia ecuaţiei ( ) va fi datã în forma implicitã de relaţia U (x, y) = C Metodã de rezolvare - se identificã în ecuaţie P(x, y) şi Q(x, y) şi se verificã egalitatea ( ) (x y) x Q x y y P , , ∂ ∂ = ∂ ∂ - se determinã funcţia U - se scrie soluţia ecuaţiei sub formã implicitã U(x, y) = C Dacã este posibil, din aceastã egalitate se aflã y în funcţie de x şi se obţine forma explicitã a soluţiei Exemplu : Sã se determine soluţia generalã a ecuaţiei (x + y + )dx + (x − y + )dy = In acest caz P(x, y) = x + y + şi Q(x, y) = x − y + şi ( , ) ( , ) = ∂ = ∂ ∂ ∂ x y x Q x y y P ( ) ( ) ( ) ( ) ( ) , , , y x xy x U x y P t dt Q x t dt t dt x t dt x y x y = + + − = ∫ + ∫ =∫ + + ∫ − + = Soluţia generalã a ecuaţiei este datã sub formã implicitã de relaţia x y + x + y + xy − = C Aceastã ecuaţie nu poate fi rezolvatã analitic în raport cu necunoscuta y , deci nu se poate preciza forma explicitã a soluţiei Ecua􀄠ii reductibile la ecua􀄠ii fundamentale Tehnica generalã de rezolvare a acestui tip de ecuaţii este urmǎtoarea : - se reduce ecuaţia la o ecuaţie fundamentalã - se rezolvã ecuaţia fundamentalã - se scrie soluţia ecuaţiei iniţiale folosind soluţia celei fundamentale Ecuatii omogene Forma generalã a unei ecuaţii omogene este y ' = f ( y / x) Prin schimbarea de variabilã ( ) ( ) y x z x x = se obţine o ecuaţie cu variabile separabile Exemple : Sǎ se determine soluţia generalǎ a ecuaţiei xy '− y = x + y Pentru x ≠ ecuaţia se scrie ' y y y x x = + + ⎛ ⎞ ⎜ ⎟ ⎝ ⎠ Cu schimbarea de variabilã ( ) ( ) y x z x x = , adicã y (x) = x ⋅ z (x) , ecuaţia devine z (x) + xz '(x) = z (x) + + z (x) , care este o ecuaţie cu variabile separabile, anume ' z x z = + Integrala generalã a ecuaţiei este dz dx z x = +∫∫ , adicã ( ) ( ) ln z + + z = ln | x | +C = ln C | x | Rezultã z + + z = Cx , adicǎ ( ) | | C x C x z x = − Soluţia generalã a ecuaţiei este ( ) : { } , | | C x y R R y x x C x − − → = Ecua􀄠ii reductibile la ecuatii omogene sau cu variabile separabile Ecuaţiile având forma generalã ' ' ' ' ax by c y f a x b y c ⎛ + + ⎞ = ⎜ ⎟ ⎝ + + ⎠ pot fi reduse la ecuaţii omogene sau cu variabile separabile astfel : - dacã a / a' ≠ b / b' se rezolvã sistemul de ecuaţii ' ' ' ax by c a x b y c + + = ⎧⎨ ⎩ + + = care are soluţia ( ) x , y Prin schimbarea de variabile x = u + x , y = v + y se obţine o ecuaţie omogenã cu variabila independentã u şi funcţia necunoscutã v - dacã a / a' = b / b' se foloseşte substituţia z = ax + by şi ecuaţia se transformã într-o ecuaţie cu variabile separabile Exemple: ( x + y − ) − (x − y − ) y ' = Ecuaţia se scrie sub forma ' x y y x y + − = − − deoarece y = x − nu este soluţie a ecuaţiei Sistemul x y x y + − = ⎧⎨ ⎩ − − = are solutia unicã x y = ⎧⎨ ⎩ = − Se face substituţia x u y v = + ⎧⎨ ⎩ = − şi se obţine ecuaţia omogenã ( u + v) + (v − u)v ' = cu funcţia necunoscutã v Notând v z u = , adicã v = zu ecuaţia se reduce la ecuaţia cu variabile separabile ' z z z u z + + = ⋅ − Integrala generalã a acestei ecuaţii este z dz du z z u − = + + ∫ ∫ Calculând cele douã integrale obţinem ( ) ( ) ln ln z z arctg z u C + + − + = + Tinând cont cã v y z u x + = = − se obţine soluţia generalã sub formã implicitã (( ) ( )( ) ( ) ) ln x y y x y x arctg x + − + + − + + − − = − Forma explicitã a soluţiei nu se poate determina ( x + y + ) − ( x + y − ) y ' = Ecuaţia se scrie sub forma ( ) ' x y y x y + + = + − Deoarece ' ' a b a b = = = se va folosi substituţia x + y = z Din z x y − = rezultã ' ' z y − = Ecuaţia devine ' z z z − + = − adicã ' z z z = − Aceasta este o ecuaţie cu variabile separabile care se poate scrie sub forma ' z z ⎛ − ⎞ = ⎜ ⎟ ⎝ ⎠ Integrala generalã a acestei ecuaţii conduce la relaţia ln z − z = x + C Tinând cont de expresia lui z se obţine soluţia generalã a ecuaţiei iniţiale, soluţie scrisã sub formã implicitã : ( ) ( ) x + y − ln x + y − x = C Nici în acest caz nu se poate preciza forma explicitã a soluţiei Ecua􀄠ii ce admit factor integrant Au forma generala P(x, y)dx +Q(x, y)dy = cu P Q y x ∂ ∂ ≠ ∂ ∂ dar pentru care existã funcţia μ = μ (x, y) ≠ , numitã factor integrant, astfel încât ( P) ( Q) y x μ μ ∂ ∂ = ∂ ∂ Dacã factorul integrant μ (x, y) poate fi determinat, atunci ecuaţia μ (x, y)P(x, y)dx +μ (x, y)Q(x, y) = este o ecuaţie cu diferenţialã totalã, echivalentã cu cea iniţialã Nu toate ecuaţiile au factor integrant, existã doar câteva cazuri importante dintre care menţionǎm: - dacã Q ∂P / ∂y − ∂Q/ ∂x depinde doar de x atunci existã factor integrant ce depinde doar de x si μ = μ (x) satisface ecuaţia / / ' P y Q x Q μ μ ∂ ∂ −∂ ∂ = ( ) - dacã P ∂Q/ ∂x − ∂P / ∂y depinde doar de y atunci atunci existã factor integrant ce depinde doar de y si μ = μ ( y) satisface ec / / ' Q x P y P μ μ ∂ ∂ −∂ ∂ = ( ) Pentru rezolvarea ecuaţiilor cu factor integrant se parcurg urmãtoarele etape : - se determinã factorul integrant rezolvând ecuaţiile diferenţiale ( ) sau ( ) - se scrie ecuaţia cu diferenţiale totale corespunzãtoare - se rezolvã ecuaţia cu diferenţiale totale (cu necunoscuta y = y(x) ) şi se obţine astfel soluţia ecuaţiei iniţiale Exemplu : ( x + y + y )dx + ( xy + x)dy = In acest caz P(x, y) = x + y + y şi Q(x, y) = xy + x Rezultã cã P y y ∂ = + ∂ şi Q y x ∂ = + ∂ Ecuaţia nu are diferenţialã totalã deoarece P Q y x ∂ ∂ ≠ ∂ ∂ Totuşi P Q y x Q x ∂ − ∂ ∂ ∂ = depinde numai de x , deci se poate alege un factor integrant de forma μ = μ (x) El va satisface ecuaţia ' x μ = μ ⋅ care este o ecuaţie cu variabile separabile cu soluţia μ (x) = x Din înmulţirea cu x a ecuaţiei iniţiale se obţine ecuaţia cu diferenţialã totalã ( x + x y + x y)dx + ( x y + x )dy = - Funcţia ( ) ( ) , x y U x y = ∫ t dt + ∫ x t + x dt = x + x y + x y Soluţia ecuaţiei, scrisã sub formã implicitã va fi deci x + x y + x y = C Ea este şi soluţia ecuaţiei iniţiale Ecuatii de tip Bernoulli Forma generalã este y'= P(x)⋅ y + Q(x)⋅ yα , unde α ∈R , α ≠ , α ≠ şi P,Q: I →R sunt funcţii date, continue pe I Pentru a > ecuaţia are soluţia singularǎ y : I →R , y(x) = Prin schimbarea de funcţie z = y −α se obţine o ecuaţie liniarã Dacǎ soluţia ecuaţiei liniare este g z atunci soluţia ecuaţiei iniţiale este = α − g y z Exemplu : y ' y x y x=+ In acest caz α = / Se foloseşte substituţia z = y − / = y / Rezultã y = z şi y ' = ⋅ z ⋅ z ' Ecuaţia devine z z ' z xz x ⋅ ⋅ = + , adicã z z ' z x x ⎛ − − ⎞ = ⎜ ⎟ ⎝ ⎠ Din soluţia z = rezultã soluţia singularã y = Ecuaţia liniarã ' x z z x = + are soluţia ( ) ln z x x K x = ⎛ + ⎞ ⋅ ⎜ ⎟ ⎝ ⎠ care conduce la ( ) ln y x x K x = ⎛ + ⎞ ⋅ ⎜ ⎟ ⎝ ⎠ Ecuatii de tip Ricatti Forma generala este y'+P(x)⋅ y + Q(x)⋅ y + R(x) = Aceste ecuaţii se pot rezolva numai dacã se cunoaşte mãcar o soluţie particularã a lor : - dacã se cunoaşte o soluţie ( ) y x , prin transformarea y = y + / z se obţine o ecuaţie liniarã şi neomogenã ; - dacã se cunosc douã soluţii y şi y , prin schimbare de variabilǎ y y z y y − = − se obţine o ecuaţie liniarã şi omogenã ; - dacã se cunosc trei soluţii y , y , y atunci soluţia se obţine direct din relaţia : y y y y C y y y y − − = − − Exemplu : Sã se rezolve ecuaţia ' x x y y y x x x + − − = − − − a) ştiind cã admite soluţia y = −x ; b) ştiind cã admite soluţiile ( ) y x = −x şi ( ) y x = − / x ; c) ştiind cã admite trei soluţii ( ) y x = −x , ( ) y x = − / x şi ( ) y x = x + a) Dacã se cunoaşte numai soluţia y se face schimbarea de variabilã y x z =− adicã y ' z ' x z =− − Se obţine ecuaţia ( ) x z ' x x x x x z z z − ⎛ − − ⎞ + ⎛ − ⎞ − ⎛ − ⎞ − = ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ ⎝ ⎠ ⎝ ⎠ din care, dupã efectuarea calculelor rezultã ecuaţia liniarã ' x z z x x + − = − − cu soluţia k x z x + = − Rezultã kx y x k − − = + b) Dacã se cunosc douã soluţii se face substituţia / y x z y x + = + adicã ( ) z x y x z − = − şi ( ) ( ) ( )( ) ( ) ' ' ' z x x z z x z xz y x z − − − − − − = − Introducând aceste expresii în ecuaţia diferenţialã obţinem (dupã calcule) ecuaţia liniarã ' z z x = care are soluţia z = cx Rezultã c x y cx − = − Observãm ca soluţia obţinutã coincide cu cea de la a) dacã considerãm c = − / k c) dacã se cunosc y , y , y , soluţia generalã se obţine direct din formula : y y y y k y y y y − − = − − de unde rezultã x k y kx − = − Ecuatii de tip Lagrange Forma generalã este y = x ⋅ A( y ') + B( y ') , unde A( y ') ≠ y ' Se deriveazã ecuaţia şi se noteazã y ' = p Se obţine o ecuaţie liniarã cu funcţia necunoscutã x şi variabila independentã p Aceastã ecuaţie are solutia de forma x = x( p) iar soluţia generalã a ecuaţiei Lagrange se dã în formã parametricã ( ) ( ) ( ) ( ) x x p y x p A p B p ⎧ = ⎪⎨ ⎩⎪ = + Exemplu : Sã se rezolve ecuaţia ( ) y = x ⋅ y ' − y ' Prin derivarea ecuaţiei se obţine ( ) y ' = y ' + xy ' y ''− y '' Se noteazã y ' = p şi se ajunge la ecuaţia p − p = ( px − ) p' în care p este funcţie de x Dacã se considerã x ca funcţie de p (se inverseazã aplicaţia p ) şi se ţine cont de faptul cã x ' = / p' (din formula de derivare a funcţiei inverse) se obţine ecuaţia liniarã ( ) ' x x p pp + − = − − pentru p( p − ) ≠ Rezultã x = (C + ln p)/ ( p − ) şi soluţia ecuaţiei este datã parametric prin x = (C + ln p)/ ( p − ), y = p (C + ln p)/ ( p − ) − p Pentru p = şi p = se obţin douã soluţii singulare : y = K şi y = x + L Înlocuind aceste funcţii în ecuaţia iniţialã se obţine K = , respectiv L = − Deci soluţiile particulare vor fi y = şi y = x − Ecuatii de tip Clairaut Ecuaţiile de tip Clairaut au forma generalã y = xy '+ B( y ') Notând y'= p ecuaţia devine y = xp + B(p) Prin derivarea sa se obţine ecuaţia p'⋅(x + B'(p)) = Dacã p'(x) = se obţine soluţia (generalã) y(x) = Cx + B(C) Din egalitatea x + B'(p) = se obţine soluţia singularã ( ) ( ) ( ) ' ' x B p y B p p B p ⎧ = − ⎪⎨ ⎩⎪ = − + scrisã sub formã parametricã Exemplu : ( ) y = xy '− y ' Soluţia generalã este y = Cx −C şi o soluţie particularã este datã parametric de x p y xp p = ⎧⎪⎨ ⎩⎪ = − Soluţia singularǎ scrisǎ sub formǎ explicitǎ este x x x y = x ⋅ − = Exerci􀄠ii propuse Sã se rezolve urmãtoarele ecua􀄠ii diferen􀄠iale sau probleme Cauchy: y '− y / x = R : y = Cx + x y '− y / x = x R : y = x / + C / x xy '+ y − ex = , y (a) = b R : y = ex / x − (ab − ea )/ x ( ) ⎪⎩ ⎪⎨ ⎧ = − − = − − ' y x x y y R : (x x x) x x y arcsin ( ) − + − + = y 'cos x + y = tgx, y ( ) = R : y = x / cos x , x∈[ ,π / ) xy '− y = y R : y = Cx / −C x (x − y) y − x y ' = R : y = / (ln | x | +C) ( − x ) y '+ xy = ax R : y = a + C x − xy '− y = x / R : y = x / + Kx y '− xy = x R : ( ) / x y x Ce− = − + xy '− y = ln x R : y = −ln x − /( x) + Cx ( x + xy )dx + ( x y + y )dy = R : x + x y + y = C (x + y)dx + (x + y)dy = R : x + xy + y = C xy ' = y, y ( ) = R : y = xy ' = y, y ( ) = R : y = x CAPITOLUL Aplica􀄠ii ale ecua􀄠iilor diferen􀄠iale de ordinul I Suprafa􀄠a de echilibru a unui lichid în rota􀄠ie (integrare directă) Problema: Un tub vertical cilindric cu raza r se roteşte rapid în jurul axei sale cu viteza unghiulară constantă ω Să se determine ecua􀄠ia intersec􀄠iei suprafe􀄠ei de rota􀄠ie a lichidului cu un plan care con􀄠ine axa tubului Variabila independentǎ, x, indicǎ distanţa unui punct faţǎ de axa de rotaţie (axa tubului) Funcţia necunoscutǎ, y = y(x) , unde y :[−r, r]→R reprezintǎ înǎlţimea lichidului în punctele aflate la distanţa x de axa cilindrului şi descrie intersecţia suprafeţei de rotaţie a lichidului cu un plan ce conţine axa tubului Figura Suprafa􀄠a de echilibru a unui lichid în rota􀄠ie Sub acţiunea forţei de inerţie lichidul se ridicǎ spre pereţii tubului y = y(x) A α c F F β G Asupra punctului de masa m situate la distanţa x de axǎ acţioneazǎ douǎ forţe: greutatea G = −mg ⋅ j şi forţa centrifugǎ Fc = mω x ⋅ i Deoarece viteza de rotaţie e constantǎ suprafaţa lichidului e stabilǎ şi forţa rezultantǎ F este perpendicularǎ pe planul tangent la suprafaţa lichidului, adicǎ unghiurile α şi β sunt egale Rezultǎ cǎ tgα = tgβ adicǎ ( ) x g y x = ⋅ ' ω Ceea ce reprezintǎ ecuaţie diferenţialǎ ce descrie curbǎ de rotaţie Prin integrare directǎ se obţine ( ) = ∫ = x +C g xdx g y x ω ω Deci curba de rotaţie este o parabolǎ Constanta C se determinǎ folosind faptul cǎ volumul lichidului este constant Dacǎ înainte de începerea rotaţiei înǎlţimea lichidului era h atunci conservarea volumului se scrie astfel ( ) ∫ ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ ⋅ ⋅ = = ⋅ + ⋅ + r k r r h y x dx r C kr C π π π unde g k ω = , ceea ce reprezintǎ o ecuaţie cu necunoscuta C Evolutia unei popula􀄠ii (variabile separabile) a Evolu􀄠ia unei popula􀄠i într-un mediu cu resurse nelimitate Problem􀓽: S􀓽 se descrie evolu􀄠ia unei popula􀄠ii a c􀓽rei creştere în fiecare moment este propor􀄠ional􀓽 cu valoarea sa în acel moment Variabila independentǎ este timpul t iar funcţia necunoscutǎ este y :[ ,+∞)→R Ea caracterizeazǎ mǎrimea populaţiei ( y(t ) reprezintǎ numǎrul de indivizi, densitatea unei culturi bacteriene etc la momentu de timp t) In timpul Δt populaţia se modificǎ cu Δy = k ⋅ y(t )⋅Δt Variaţia sa instantanee este Δy / Δt deci se poate scrie ( ) k y(t ) t y y t t = ⋅ Δ = Δ Δ → ' lim ceea ce reprezintǎ ecuaţia diferenţialǎ ce descrie evoluţia populaţiei Ea este o ecuaţie cu variabile separabile Soluţia singularǎ y R s :[ ,+∞)→ , y (t ) = s care nu convine problemei deoarece în aceastǎ situaţie populaţia nu existǎ Soluţia generalǎ a ecuaţiei este y(t ) = C ⋅ ekt Aceasta este faimoasa lege a lui Malthus care a prevǎut cǎ populaţia Terrei va creşte exponenţial în timp ce resursele sale cresc in progresie geometricǎ, deci ele vor deveni insuficiente pentru toti oamenii şi rǎzboaiele sunt necesare pentru reglementarea echilibrului Problem􀓽 Dac􀓽 popula􀄠ia unei 􀄠􀓽ri s-a dublat în de ani, peste cât timp se va tripla, dac􀓽 viteza de creştere este propor􀄠ional􀓽 cu m􀓽rimea sa? Tinând cont cǎ y(t ) = Cekt şi cǎ y( ) = y( ) rezultǎ Ce k = C , deci ln k = Triplarea populaţiei înseamnǎ y(t ) = y( ) ceea ce conduce la ekt = , adicǎ ln ln ln = = ⋅ ≈ k t Deci populaţia se va tripla dupa de ani Problem􀓽: Intr-o colonie de microbi care se înmul􀄠esc cu o vitez􀓽 propor􀄠ional􀓽 cu num􀓽rul lor se constat􀓽 c􀓽 num􀓽rul de microbi sa dublat în ore Ce se va întâmpla dup􀓽 ore? Din legea de evoluţie y(t ) = C ⋅ ekt şi condiţia y( ) = y( ) rezultǎ ln k = Atunci y( ) = Ce ln = C aratǎ cǎ numǎrul de microbi a crescut de ori în ore b Evolu􀄠ia unei popula􀄠ii într-un mediu cu resurse limitate Problem􀓽: S􀓽 se descrie evolu􀄠ia unei popula􀄠ii într-un mediu cu resurse limitate în care popula􀄠ia maxim􀓽 ce poate tr􀓽i are m􀓽rimea L, dac􀓽 m􀓽rimea ei ini􀄠ial􀓽 este P Folosind notaţiile din problema anterioarǎ obţinem problema Cauchy ( ) ( ) ( ( )) ( ) ⎩ ⎨ ⎧ = = ⋅ − ' y P y t ky t L y t Ecuaţia diferenţialǎ, numitǎ şi ecuaţia logisticǎ, aratǎ cǎ o populaţie va creşte (adicǎ y'(t ) > ) doar dacǎ ea nu a atins nivelul maxim permis de mediu Dacǎ mǎrimea ei este mai mare decât acest nivel ea va trebui sǎ scadǎ pânǎ la limita acceptabilǎ Ecuaţia are variabile separabile Soluţile singulare y y R s s , :[ ,∞)→ , ( ) y t = s , respectiv y (t ) L s = descriu situaţia în care populaţia nu existǎ (adicǎ P = ), respective siuaţia când populaţia iniţialǎ este la nivelul maxim admis de mediu, adicǎ P = L şi rǎmane la acest nivel Scrisǎ sub forma ( ) ( )( ( )) k y t L y t y t = − ' ecuaţia poate fi integratǎ şi conduce la ( ) ( ) Lkt C L y t y t + = ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ − ln , adicǎ ( ) ( ) CeLkt L y t y t = − Soluţia gneralǎ este deci ( ) Lkt Lkt Ce Ce y t L + = Determinarea constantei C se face folosind condiţia iniţialǎ Din ( ) C LC y P + = = rezultǎ P P C − = deci soluţia problemei este ( ) Lkt Lkt L P P e L P e y t − + ⋅ ⋅ = Din expresia analiticǎ a soluţiei se observǎ cǎ y(t ) L t = →∞ lim , deci populaţia tinde sǎ atingǎ valoarea maximǎ admisǎ Varia􀄠ia presiunii atmosferice în raport cu altitudinea (ecua􀄠ie cu variabile separabile) Problem􀓽: S􀓽 se determine presiunea atmosferic􀓽 p = p(h) a unei mase volumice de aer ρ la temperature T în raport de în􀓽l􀄠imea h m􀓽surat􀓽 de la nivelul m􀓽rii, daca presiunea la nivelul m􀓽rii este ( ) p = p Variabila independentǎ este altitudinea h iar funcţia necunoscutǎ este presiunea atmosfericǎ p = p(h) Ecuaţia fundamentalǎ a hidrostaticii, scrisǎ pentru acest caz, este p'+ρg = Dacǎ vom considera aerul la un gaz perfect obţinem relaţia p ⋅ v = n ⋅ R ⋅T , unde n este numǎrul de moli de aer ce ocupǎ volumul v, R este o constantǎ Rezultǎ cǎ R T M p n R T m p v m ⋅ = ⋅ ⋅ ⋅ ρ = = ⋅ , unde M este masa molarǎ a aerului In acest caz ecuaţia ce descrie presiunea aerului devine ' = ⋅ + ⋅ p R T M g p Problema Cauchy ataşatǎ este ( ) ( ) ( ) ( ) ' ⎪⎩ ⎪⎨ ⎧ = = ⋅ + ⋅ p p p h R T h M g p h Cazul I Dacǎ temperatura este constantǎ, ecuaţia se scrie p'= −kp , unde R T M g k ⋅ = ⋅ In acest caz ecuaţia are variabile separabile Soluţia singularǎ p(h) = nu convine problemei şi din soluţia generalǎ p(h) = C ⋅ e−kh deducem, impunând condiţia iniţialǎ cǎ soluţia problemei Cauchy este p(h) = p ⋅ e−kh Rezultǎ deci cǎ presiunea scade exponenţial în raport cu altitudinea Cazul II: Dacǎ temperatura nu este constantǎ dar aerul respectǎ legea transformǎrilor adiabatice p ⋅ vγ = K , cu γ ≠ dat, atunci /γ ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ = p K v , din relaţia p ⋅ v = n ⋅ R ⋅T se obţine n K p p K n p v n p R T /γ /γ ⋅ − /γ = ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ ⋅ = = şi problema Cauchy devine ( ) ( ) ( ) ( ) ⎪⎩ ⎪⎨ ⎧ = = − ⋅ = − ⋅ / / / ' p p p h k p h K m g p h notatie γ γ γ Integrând ecuaţia ( ) ( ) p− /γ h ⋅ p' h = −k obţinem ( ) ( ) C h k h p + − ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ − = − / γ γ , adicǎ ( ) ( ) − ⎥⎦ ⎤ ⎢⎣ ⎡ + − ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ = − γ γ γ p h k h C Din condiţia iniţialǎ rezultǎ (γ ) γ γ γ / − − C = p , deci ( ) ( ) / − − ⎥⎦ ⎤ ⎢⎣ ⎡ ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ − = − γ γ γ γ γ γ p h p k h ceea ce aratǎ cǎ presiunea scade atunci când creşte altitudinea şi altitudinea maximǎ la care existǎ atmosferǎ este ( ) ( ) / max − = ⋅ − γ γ γ γ p h care se obţine din relaţia ( ) max p h = Căderea liberă (cu variabile separabile) Problem􀓽: S􀓽 se determine viteza unui corp în mişcare vertical􀓽 sub ac􀄠iunea greut􀓽􀄠ii sale şi a rezisten􀄠ei aerului, dac􀓽 viteza de pornire este v Variabila independentǎ este timpul si funcţia necunoscutǎ este viteza v = v(t ) Condiţia iniţialǎ este ( ) v = v Rezistenţa aerului este R=R(t) şi acceleraţia corpului este a(t ) = v'(t ) Legea fundamentalǎ a dinamicii se scrie sub forma m⋅ g − R(t ) = mv' Cazul I: Dacǎ rezistenţa aerului este proporţionalǎ cu viteza corpului, adicǎ R(t ) = kv(t ), ecuaţia devine mg − kv(t) = mv'(t ) Notând k m/ k = obţinem problema Cauchy () () ( ) ' ⎩ ⎨ ⎧ = = − v v v t g k v t Ecuaţia diferenţialǎ, cu variabile separabile, are soluţia singularǎ ( ) v t = g / k Pentru determinarea soluţiei generale se scrie ecuaţia sub forma ( ) ( ) ' = g − k v t v t , deci ( ) ( ) t C g k v t v t = + ∫ − ' , adicǎ g k v(t ) t C k − ln( − ) = + De aici rezultǎ ( ) ( C e k t ) k g v t = − − Din condiţia iniţialǎ ( ) C k g v = − rezultǎ g k v C = − In aceastǎ situaţie viteza creşte odatǎ cu trecerea timpului şi tinde sǎ atingǎ valoarea maximǎ posibilǎ, max v = g / k Cazul II Dacǎ rezistenţa aerului este proporţionalǎ cu pǎtratul vitezei (ceea ca se întâmplǎ când vitezele de mişcare sunt mari) atunci R(t ) = kv (t ) şi problema Cauchy devine ( ( )) ( ) ⎩ ⎨ ⎧ = = − ' v v v g α v t unde α = k /m Ecuaţia diferenţialǎ cu variabile separabile are soluţia singularǎ v(t ) = /α care este o funcţie constantǎ şi corespunde cazului când v = /α Soluţia generalǎ, obţinutǎ prin integrarea ecuaţiei ( ) ( ) g v t v t = − ' α , este datǎ de egalitatea ( ) ( ) gt C v t v t + = ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ − + α α α ln Rezultǎ ( ) ( ) ( ) + − = ⋅⋅ ⋅ + ⋅⋅ ⋅ + g t C g t C e e v t α α α Valoarea lui C se determinǎ din condiţiile iniţiale Din analiza expresiei lui v se poate deduce cǎ viteza creşte pe mǎsurǎ ce timpul trece şi tinde sǎ atingǎ valoarea maximǎ admisǎ v / m/ k max = α = Descărcarea unui condensator într-o resisten􀄠ă (liniară) Fenomenul este important pentru cã apare în circuitele folosite la transmisiunile radio, de televiziune, la radare etc Problema : Se considerã un circuit electric format dintr-un condensator cu capacitatea C şi o rezisten􀄠ã R Se cere intensitatea curentului, i (t ) şi diferen􀄠a de poten􀄠ial v(t ) la bornele condensatorului în func􀄠ie de momentul t la care se face mãsurarea dacã sarcina ini􀄠ialã este Q Rezolvare : Considerãm funcţiile i,q,v :[ ,+∞)→R cate indicǎ intensitatea, sarcina electricǎ şi diferenţa de potenţial la bornele condensatorului în funcţie de timpul t Intre ele existã relaţia q(t ) = Cv(t ) Intensitatea curentului electric la descãrcare este i (t ) = −q'(t ) şi la bornele rezistenţei ea satisface relaţia v i R = (legea lui Ohm) Relaţiile anterioare aratã cã acest circuit este caracterizat de problema Cauchy ( ) ( ) ( ) ' q t q t C R q Q ⎧ ⎪ = − ⋅ ⎨⎪ ⎩ = (CR) în care q = q(t ) este funcţia necunoscutã iar C şi R sunt constante date în problemã Ecuaţia diferenţialã poate fi consideratã ca o ecuaţie cu variabile separabile sau ca o ecuaţie liniarã şi omogenã de ordinul I Soluţia problemei Cauchy este ( ) CR q t Qe − = Intensitatea curentului (la descãrcare) este ( ) ( ) /( ) ' t CR t CR Q i t q t e I e CR − = − = = − iar ( ) ( ) t CR Q v t R i t e C − = ⋅ = Momentul începând de la care descãrcarea este practic terminatã se considerã a fi τ pentru care ( ) I i τ = Se obţine eCR −τ = adicã τ = ⋅C ⋅ R ⋅ ln ≈ , ⋅C ⋅ R Incãrcarea unui condensator printr-o rezisten􀄠ã în prezen􀄠a unei surse de curent continuu (liniară) Problemã : Se considerã un circuit alcãtuit dintr-un condensator cu capacitatea C , o rezisten􀄠ã R şi o sursã de curent continuu având for􀄠a electromotoare constantã E Se cere sã se determine intensitatea curentului şi diferen􀄠a de poten􀄠ial la bornele condensatorului în func􀄠ie de momentul la care se face mãsurarea Rezolvare : Se considerã funcţiile i,q,v :[ ,+∞)→R care reprezintã intensitatea curentului, sarcina condensatorului şi diferenţa de potenţial la bornele acestuia în funcţie de timp La încãrcarea condensatorului i (t ) = q'(t ) iar ( ) q(t ) v t C = Legea lui Kirchoff aratã cã ( ) ( ) E v t i t R − = , deci problema Cauchy ce caracterizeazã circuitul este ( ) ( ) ( ) ' q t R q t E C q ⎧ ⋅ + = ⎪⎨⎪ ⎩ = Soluţia generalã a ecuaţiei (liniarã şi neomogenã de ordinul I) este ( ) t CR q t CE Ke − = + iar soluţia problemei Cauchy este ( ) t CR q t C E e − ⎛ ⎞ = ⋅ ⋅ ⎜⎜ − ⎟⎟ ⎝ ⎠ Rezultã imediat ( ) ( ) t CR q t v t E e C ⎛ − ⎞ = = ⎜⎜ − ⎟⎟ ⎝ ⎠ şi ( ) '( ) t CR E i t q t e R − = = Momentul în care condensatorul este practic încãrcat este cel la care diferenţa de potenţial este ( ) v τ = E Din E E e CR τ ⎛ − ⎞ = ⎜⎜ − ⎟⎟ ⎝ ⎠ rezultã τ = C ⋅ R ⋅ ln ≈ , ⋅C ⋅ R Transformarea energiei electrice in căldură (liniară) Problem􀓽: S􀓽 se descrie modificarea temperaturii unui corp în raport cu mediul ambiant atunci când pentru înc􀓽lzire se foloseşte energia electric􀓽 Variabila independentǎ este timpul t (mǎsurat în secunde) iar funcţia necunoscutǎ este temperatura corpului θ =θ (t ) presupusǎ aceeaşi în toate punctele corpului Se considerǎ W puterea electricǎ mǎsuratǎ în waţi, m masa corpului mǎsuratǎ în grame, c cǎldura masicǎ (cantitatea de cǎldurǎ necesarǎ pentru a ridica temperatura unui gram de material cu un grad), S suprafaţa de rǎcire, α coeficientul de împrǎştiere (cantitatea de cǎdurǎ împrǎştiatǎ de cm într-o secundǎ pentru ridicarea temperaturii cu un grad) Pentru obţinerea ecuaţiei se foloseşte principiul conservǎrii energiei, considerând cǎ în durata de timp Δt energia electricǎ este utilizatǎ pentru a ridica temperatura corpului (cu masa m grame) cu Δθ în timp ce o parte din cǎldurǎeste împrǎştiatǎ: Relaţia m c S t W t = ⋅ ⋅Δ + ⋅ ⋅ ⋅Δ ⋅Δ θ α θ conduce la ecuaţia diferenţialǎ m c W m c S ⋅ ⋅ ⋅ = ⋅ + ⋅ ' θ θ α Se noteazǎ m c S B ⋅ =α ⋅ şi m c W D ⋅ ⋅ = şi ecuaţia devine Este o ecuaţie liniarǎ de ordinul I, dar poate fi socotitǎ şi ecuaţie cu variabile separabile Ecuaţia omogenǎθ '(t )+ B⋅θ (t ) = are soluţia generalǎ θ (t ) = C ⋅ e−Bt Aplicând metoda variaţiei constantei considerǎm θ (t ) = C(t )⋅ e−Bt şi ecuaţia devine C'(t ) = D⋅ eBt Rezultǎ ( ) e K B D C t = Bt + Rezultǎ cǎ soluţia generalǎ a ecuaţiei este ( ) / Bt K e Bt S W t D B K e− + ⋅ − ⋅ ⋅ = + ⋅ = α θ Dacǎ temperatura iniţialǎ θ este precizatǎ atunci S W K ⋅ ⋅ = − α θ In aceastǎ situaţie temperatura corpului va fi ( ) ( ) e Bt e Bt S W t − − + − ⋅ ⋅ = θ α θ Temperatura maximǎ a corpului va fi ( ) S W t t ⋅ ⋅ = = →∞ α θ θ lim max , indiferent de temperatura de pornire Dacǎ nu se foloseşte electricitatea pentru încǎlzire atunci W= şi temperatura corpului este (t ) = e−Bt θ θ Ea descreşte exponenţial la Formula fundamentală a curentului alternativ (liniară) Problemã : Se considerã un circuit în care ac􀄠ioneazã o for􀄠ã electromotoare datoratã unei varia􀄠ii de flux şi con􀄠inând o rezisten􀄠􀓽 R şi o bobinã cu inductan􀄠a proprie L montate în serie Sã se determine intensitatea curentului electric din circuit Rezolvare : Pentru a obţine o variaţie a fluxului electric Φ(t ) se considerã un cadru cu n spire de arie S mişcându-se într-un câmp magnetic cu inducţia B , cadru închis printr-un circuit exterior Acest cadru se roteşte uniform cu viteza unghularã ω Fluxul captat Φ(t ) se descompune în douã pãrţi : - Fluxul ( ) Φ t = n ⋅ B ⋅ S ⋅ sinωt provenind de la polul nord al câmpului magnetic - Fluxul ( ) ( ) Φ t = L ⋅ i t generat de cadrul parcurs de curentul electric Din relaţia (datã în problemã) E(t ) = −Φ'(t ) şi ţinând cont de faptul cã ( ) E(t ) i t R = se obţine ecuaţia ( ) cos '( ) n B S L i t t i t R R ω ω ⋅ ⋅ ⋅ = − − ⋅ Ea se scrie sub forma L ⋅ i '(t ) + R ⋅ i (t ) = −n ⋅ B ⋅ S ⋅ω ⋅ cosωt = E cosωt Din rezolvarea ecuaţiei omogene se obţine i(t ) = C ⋅ e−Rt / L Aplicând metoda variaţiei constantei obţinem t L E Ce L R e L R C'e Rt / L C RT / L + Rt / L = cosω − + − − − , deci ( ) e t L E C t t L R ' = ⋅ cosω Rezultǎ ( ) ∫ ( − ) + + = = R t L t e k R L E e tdt L E C t Rt L Rt / L / cos sinω ω cosω ω ω Soluţia generalǎ a ecuaţiei este ( ) (R t L t ) ke Rt L R L E i t / sin − cos + − + = ω ω ω ω Aceasta are o componentã periodicã dominantã (anume Rcosωt + Lω sinωt ) şi o componentã neglijabilã (anume R t Re L − ) atunci când timpul t este mare Din acest motiv curentul obţinut se numeşte curent alternativ Oglinda parabolică (ecua􀄠ie omogenă) Problem􀓽: S􀓽 se determine forma unei oglinzi care reflect􀓽 o raz􀓽 luminoas􀓽 ce porneşte din origine pe o direc􀄠ie paralel􀓽 cu Ox Graficul funcţiei necunoscute y=y(x) reprezintǎ oglinda în care se reflectǎ lumina Figura Oglinda parabolicǎ Raza OM este reflectatǎ de oglindǎ pe direcţia MA Considerând tangenta MB la graficul funcţiei y=y(x), obţinem ∧ ∧ ∧ BMA = CMO = MCO deoarece unghiul de incidenţǎ este egal cu unghiul de reflexie Atunci ∧ ∧ MOP = BMA, deci ( ) ( ) ( ( )) ' ' ( ) y x y x tg BMA tg BMA tg BMA x y x tg MOP − = ⎟⎠ ⎞ ⎜⎝ − ⎛ ⎟⎠⎞ ⎜⎝ ⎛ = = = ⎟⎠ ⎞ ⎜⎝ ⎛ ∧ ∧ ∧ ∧ şi se obţine ecuaţia omogenǎ ( ) x y y y = − ' ' In mod obişnuit pentru rezolvarea acestei ecuaţii se foloseşte substituţia y / x = u , dar în acest caz calculele sunt complicate şi este mai util un artificiu de calcul datorat lui Lagrange: se noteazǎ y'= m şi se deriveazǎ ecuaţia m m y x − = ⋅ M B O P N R y = y(x) C Se obţine o nouǎ ecuaţie m m m m x m m y ⋅ = − + + − = ' ( ) ' din care rezultǎ cǎ ( ) x m m m ' − = , adicǎ o ecuaţie cu variabile separabile Soluţile singulare m(x) = şi m(x) = ± nu convin problemei pentru cǎ în acest caz se anuleazǎ numitorul din ecuaţia iniţialǎ Separând variabilele obţinem m m x m ( ) ' = − care, prin integrare conduce la (kx) m m ln ln = ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ − , adicǎ kx m m = − Introducând expresia km m x − = în ecuaţia m m y x − = ⋅ se obţine ' k m k y y ⋅ − = ⋅ = − , adicǎ ky m = − Din ( ( )) ( ( )) / / k ky ky km m x − − − = − = rezultǎ ⎟⎠ ⎞ ⎜⎝ = ⎛ − x k k y ceea ce aratǎ cǎ oglinda are forma unei parabole (se numeşte oglindǎ parabolicǎ) al cǎrei focar este originea axelor de coordonate In fapt ea reprezintǎ interiorul unui paraboloid de rotaţie Exerci􀄠ii propuse Se ştie cǎ viteza de rǎcire a unui corp este proporţionalǎ cu diferenţa de temperaturǎ între corp si mediul înconjurǎtor a) Sǎ se scrie relaţia care existǎ între temperatura T şi timpul t, dacǎ un corp încǎlzit la T grade este plasat într-o camerǎ a cǎrei temperaturǎ constantǎ este de a grade b) In cât timp va scadea la grade temperature unui corp încǎlzit la grade dacǎ temperatura camerei în care e introdus este grade şi în primele min corpul se rǎceşte la grade? R: a) T'(t ) = k(T(t )− a) are soluţia T(t ) = (T − a)e−kt b) , T = a = ; din T( ) = rezultǎ ( ) / e−k = / ; timpul necesar este min Frânarea unui disc care se roteşte uniform într-un lichid este proporţionalǎ cu viteza unghiularǎ a) Scrieţi şi rezolvaţi ecuaţia care descrie modificǎrile vitezei unghiulare în timpul rotaţiei b) Care este viteza dicului dupa min dacǎ viteza sa iniţialǎ a fost de turaţii/min şi dupǎ un minut viteza sa a scǎzut la turaţii /min R: a) ω'(t ) = −kω(t ) , k> , are soluţia ω(t ) = C ⋅ e−kt ; semnul “-“ ce apare în ecuaţie aratǎ cǎ miscarea este încetinitǎ b) ( ) t t ⎟⎠ ⎞ ⎜⎝ = ⋅⎛ ω deci ω( ) = grade Viteza de dezintegrare a radiului este proporţionalǎ cu cantitatea sa Se ştie cǎ dupǎ ani nu rǎmâne decât o jumǎtate din cantitatea de radiu iniţialǎ (perioada de înjumǎtǎţire este ani) Sǎ se determine procentajul de radiu dezintegrat dupǎ ani R: Se considerǎ Q(t) cantitatea de radiu dezintegratǎ sin cantitatea iniţialǎ C Ecuaţia liniarǎ '( ) ( ( )) Q t = k ⋅ C −Q t şi condiţia iniţialǎ Q( ) = conduc la ( ) ( ) Q t = C − e−kt Din condiţia Q( ) = ⋅Q( ) rezultǎ ( ) / e−k = / Procentajul cerut este ⋅ (Q( )/Q( ))% = ( − / )⋅ % ≈ % Viteza de scurgere a apei printr-o deschidere aflatǎ pe verticalǎ la distanţa h de suprafaţa apei este datǎ de formula v = c gh unde c ≈ şi g este acceleraţia gravitaţionalǎ In cât timp se va goli un rezervor cubic cu latura de m printr-o deschidere practicatǎ la baza sa R: Se considerǎ h = h(t ) înǎlţimea lichidului scurs pânǎ la momentul t Ecuaţia diferenţialǎ h'(t ) = c g( − h(t)) şi condiţia iniţialǎ h( ) = caracterizeazǎ fenomenul de curgere Soluţia problemei Cauchy este ( ) ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ = − t c g h t c g Rezervorul se va goli atunci când ( )= g h t , adicǎ = ≈ c g tg unitǎţi de timp Cantitatea de luminǎ absorbitǎ de un rezervor de apǎ este proporţionalǎ cu cantitatea de luminǎ incidentǎ şi cu adâncimea rezervorului Se ştie cǎ într-un rezervor cu adâncimea de m apa absoarbe jumǎtate din catitatea de luminǎ iniţialǎ Care este cantitatea de luminǎ care poate ajunge la o adâncime de m? R: Cantitatea de luminǎ absorbitǎ, Q = Q(h) depinde de adâncimea la care se aflǎ observatorul Ecuaţia diferenţialǎ ce caracterizeazǎ fenmenul de absorbţie este '( ) ( ( )) Q h = k C −Q h şi condiţia iniţialǎ impusǎ este Q( ) = Q( )/ Soluţia problemei Cauchy asociatǎ este ( ) ( ( ) / ) / h Q h = Q ⋅ − Cantitatea de luminǎ absorbitǎ pânǎ la m este ( ) ( ( / ) ) Q = Q − iar cantitatea de luminǎ ce pǎtrunde la acel nivel este ( ) Q ⋅ / Rezistenţa aerului exercitatǎ asupra unui corp lansat cu o paraşutǎ este proporţionalǎ cu pǎtratul vitezei de mişcare Sǎ se determine viteza limitǎ pe care o poate atinge un corp R: Ecuaţia diferenţalǎ este mv'= mg − kv şi are soluţia ( ) / ( ) / ( ) + = ⋅ − ⋅⋅ ⋅ + ⋅⋅ ⋅ + k m g t C k m g t C e e k m v t Valoarea lui C se determinǎ din condiţiile iniţiale Din analiza expresiei lui v se poate deduce cǎ viteza creşte pe mǎsurǎ ce timpul trece şi tinde sǎ atingǎ valoarea maximǎ admisǎ v / m/ k max = α = Baza unui rezervor de l este acoperitǎ cu un amestec de sare şi substanţǎ insolubilǎ Presupunând cǎ viteza de dizolvarea sǎrii este proporţionalǎ cu diferenţa între concentraţia instantanee şi concentraţia soluţiei saturate ( kg sare pentru l apǎ) si cantitatea de apǎ purǎ dizolvǎ / kg sare pe minut, sǎ se determine cantitatea de sare din soluţie dupǎ o orǎ R: se noteazǎ x(t ) cantitatea de sare din soluţie la momentul t (sare dizolvatǎ în apǎ) Ecuaţia diferenţialǎ este ( ) ⎟⎠ ⎞ ⎜⎝ = ⎛ − ' x x t k Soluţia problemei Cauchy care foloseşte condiţia iniţialǎ x( ) = este ( ) ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ = − − kt x t e Dupǎ o orǎ cantitatea de sare din soluţie va fi ( ) ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ = − − k x e CAPITOLUL Ecua􀄠ii diferen􀄠iale de ordin superior Ecua􀄠ii liniare O problemã inportantã este rezolvarea ecuaţiilor diferenţiale de ordin mai mare ca Sunt puţine ecuaţiile pentru care se poate preciza forma analiticã a soluţiei Cel mai frecvent utilizate sunt ecuaţiile liniare Forma generalã a ecuaţiei liniare de ordin n este ( ) ( ) ( ) ( ) ( ) n n n y a x y a x y f x + − + + = ( ) Ecuaţia liniarã omogenã asociatã ecuaţiei ( ) este ( ) ( ) ( ) ( ) n n n y a x y a x y + − + + = ( ) În legãturã cu ecuaţiile liniare şi omogene se poate demonstra urmatorul rezultat important Teorema a) Dacã p∈N şi p y , y , , y sunt soluţii ale ecuaţiei ( ) iar C C C R p , , , ∈ atunci p p y = C y +C y + +C y este soluţie a ecuaţiei ( ) b) Mulţimea soluţiilor ecuaţiei ( ) formeazã un spaţiu vectorial de dimensiune n c) Dacã ecuaţia ( ) admite soluţia complexǎ y = u + i ⋅ v atunci funcţiile reale u şi v sunt soluţii ale ecuaţiei ( ) Observa􀄠ie : pentru determinarea soluţiei generale e ecuaţiei omogene trebuie determinate n soluţii liniar independente Teorema Soluţiile y y y I R n , , , : → ale ecuaţiei ( ) sunt liniar independenta dacã şi numai dacã existã x ∈I astfel încât determinantul ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) y (x) y (x) y (x) y x y x y x y x y x y x W y y y x n n n n ' ' ' ' ' ' , , , = (numit wronskianul sistemului) sã fie nenul în x Observa􀄠ii : ) Teorema Abel-Ostrogradski-Liouville aratã cã, dacã I este un interval ce conţine pe x şi ( , , , )( ) W y y y x ≠ n , atunci ( , , , )( ) W y y y x ≠ n pentru orice x∈I ) Dacã n y , y , , y sunt soluţii liniar independente ale ecuaţiei ( ) şi C C C R n , , , ∈ atunci soluţia generalã a ecuaţiei ( ) este n n y = C y +C y + +C y ( ) În cazul sistemelor cu coeficienţi constanţi determinarea soluţiilor liniar independente se face cu ajutorul ecuaţiei caracteristice, dar reprezintã o problemã complicatã în cazul sistemelor cu coeficienţi variabili Pentru sistemele neomogene sa poate arãta cã : Teorema : Solu􀄠ia generalã a ecua􀄠iei ( ) este suma dintre solu􀄠ia generalã a ecua􀄠iei omogene ataşatã, ( ), şi o solu􀄠ie particularã a ecua􀄠iei ( ) Metoda de rezolvare a ecuaţiilor liniare are trei paşi: - se rezolvã ecuaţia omogenã şi se obţine soluţia G y - se determinã o solutie P y a ecuaţiei neomogene - se scrie soluţia generalã a ecuaţiei neomogene G P y = y + y Ecua􀄠ii liniare cu coeficien􀄠i constan􀄠i Pentru rezolvarea ecuaţiilor cu coeficienţi constanţi este cunoscut algoritmul de rezolvare a ecuaţiei omogene Dacă funcţia f (x)are anumite forme particulare atunci există reguli şi pentru determinarea unei soluţii particulare A ) Rezolvarea ecua􀄠iei omogene Forma generalã a unei ecuaţii cu coeficienţi constanţi este ( ) ( ) ' n n n n y ay a y ay − + + + − + = ( ) Problema rezolvãrii ecuaţiei ( ) se reduce deci la determinarea unul sistem fundamental de soluţii În cele ce urmeazã prezentãm principala metodã de rezolvare pentru ecuaţiile liniare O soluţie a ecuaţiei se cautã sub forma ( ) x y x e λ = , prin analogie cu cazul n = Prin înlocuire în ecuaţia ( ) se obţine, dupã simplificarea cu eλ x , ecuaţia caracteristicã n n n n λ a λ − a λ a + + + − + = ( ) Teorema : Fie , , , n λ λ λ solu􀄠iile ecua􀄠iei ( ) a) - dacã , , n λ λ λ sunt reale şi distincte ale ecua􀄠iei ( ), atunci , , , n x x x e e eλ λ λ sunt solu􀄠ii liniar independente ale ecua􀄠iei ( ) b) - dacã λ este rãdãcinã realã cu ordinul de multiplicitate p pentru ecua􀄠ia ( ), atunci e x λ , xe x λ , …, x p e x − λ sunt p solu􀄠ii liniar independenteale ecua􀄠iei ( ) c) -dacã λ = a + ib este rãdãcinã complexã de ordinul p a ecua􀄠iei ( ) atunci ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) cos , sin cos , sin cos , sin cos , sin ax ax ax ax ax ax n ax n ax e bx e bx xe bx x e bx xe bx xe bx x − e bx x − e bx sunt solu􀄠ii liniar independente ale ecua􀄠iei ( ) Sistemul fundamental de soluţii se obţine prin însumarea soluţiilor liniar independente corespunzãtoare tuturor rãdãcinilor ecuaţiei ( ), iar soluţia generalã a ecuaţiei ( ) se obţine folosind formula ( ) Exemple : Sã se determine soluţia generalã a urmãtoarelor ecuaţii : y '''− y ''+ y '− y = Ecuaţia caracteristicã este λ − λ + λ − = şi are soluţiile λ = , λ = , λ = , Se aplicã a) din Teorema şi se obţine sistemul fundamental de (trei) soluţii format din y = ex , y = e x , y = e x Soluţia generalã este ( ) y x = C ex + C e x + C e x y '''+ y ''+ y '+ y = Ecuaţia caracteristicã este λ + λ + λ + = şi are soluţiile λ =λ =λ = Se aplicã b) din Teoremã pentru p = Sistemul fundamental de (trei) soluţii este format din y = ex , y = xex , y = x ex , iar soluţia generalã este ( ) y x = C ex + C xex + C x ex y ''+ y '+ y = Ecuaţia caracteristicã este λ + λ + = şi are soluţiile λ = − + i şi λ = − − i deci se aplicã c) din Teoremã pentru a = − , b = , p = Sistemul fundamental de (douã) soluţii este format din soluţiile y = e− x cos x şi y = e− x sin x iar soluţia generalã este ( ) y x = C e− x cos x + C e− x sin x yIV − y '''+ y ''− y '+ y = Ecuaţia caracteristicã este λ − λ + λ − λ + = cu soluţiile λ =λ = , λ = i şi λ = −i Soluţia generalã este ( ) y x = C e x + C xe x + C cos x + C sin x A ) Determinarea unei solu􀄠ii particulare a ecua􀄠iei neomogene Forma generalã a ecuaţiei neomogene cu coeficienţi constanţi este ( ) ( ) ( ) ' n n n n y ay a y ay f x − + + + − + = Nu existã metode generale de determinare a unei soluţii particulare dar, în unele cazuri simple, se pot folosi rezultatele urmãtoare : Dacã f (x) = P(x) este un polinom de grad k atunci soluţia particularã este un polinom de acelaşi grad, cu coeficienţi necunoscuţi care se vor determina prin înlocuirea în ecuaţie a) Dacã f (x) = eaxP(x) unde P(x) este un polinom de grad k existã douã situaţii - Dacã a nu este rãdãcinã a ecuaţiei caracteristice soluţia particularã se cautã sub forma ax ( ) Py = e Q x , unde Q este un polinom de grad k cu coeficienţi necunoscuţi - Dacã a exte rãdãcinã de ordin r a ecuaţiei caracteristice atunci soluţia particularã se cautã sub forma r ax ( ) Py = x e Q x , unde Q este un polinom de grad k cu coeficienţi necunoscuţi b) Dacã f (x) = eax (P(x)cos (bx) +Q(x)sin(bx)) atunci existã de asemeni douã situaţii - Dacã z = a + bi nu este soluţie a ecuaţiei caracteristice soluţia particularã se cautã sub forma ( ) ax ( ( )cos ( ) ( )sin ( )) p yx = e S x bx + T x bx , unde R(x) şi S ( x) sunt polinoame cu coeficienţi necunoscuţi avand drept grad cel mai mare dintre gradele lui P şi Q - Dacã z = a + bi este soluţie de ordin r a ecuaţiei caracteristice soluţia particularã se cautã sub forma ( ) r ax ( ( )cos ( ) ( )sin ( )) p y x = x e S x bx + T x bx , unde T(x) şi S (x) sunt polinoame cu coeficienţi necunoscuţi avand drept grad cel mai mare dintre gradele lui P şi Q În unele situaţii, pentru determinarea soluţiilor particulare, se poate aplica principiul superpoziţiei : Dacã P y şi P y sunt soluţii ale ecuaţiilor ( ) ( ) ( ) ( ) ( ) n n n y a x y a x y f x + − + + = , respectiv y( ) a (x)y( ) a (x)y g(x) n n + n− + + = , atunci P! P y + y este soluţie a ecuaţiei y( ) a (x)y( ) a (x)y f (x) g(x) n n + n− + + = + Exemple : Sã se determine câte o solu􀄠ie particularã pentru urmãtoarelor ecua􀄠ii : y ''− y '+ y = x Funcţia f (x) este un polinom de gradul I, deci soluţia particularã se cautã sub forma ( ) P y x = ax + b Atunci y '(x) = a şi y ''(x) = Introducând în ecuaţie obţinem − a + ax + b = x şi din identificarea coeficienţilor rezultã a = şi, adicã b = Soluţia particularã este ( ) P y x = x + Soluţia generalã a ecuaţiei este ( ) y x = C ex + C xex + x + y ''− y '+ y = xex Deoarece f (x) = xex = e ⋅xP(x) se încadreazã la b) pentru a = , P(x) = x şi a = este rãdãcinã dublã a ecuaţiei caracteristice , soluţia particularã se va cãuta sub forma ( ) x ( ) P y x = x e Ax + B Atunci y '(x) = ex (Ax + Bx + Ax + x) şi y ''(x) = ex (Ax + Bx + Ax + Bx + Ax + B) Introducând în ecuaţie obţinem ex ( Ax + B) = xex Din identificarea coeficienţilor se obţine A = / şi B = Soluţia particularã este deci x / Py = x e Soluţia generalã a ecuaţiei este ( ) y x = C ex + C xex + x ex / y ''+ y = xsin x Funcţia f (x) = xsin x = e x ( ⋅ cos x + xsin x) se încadreazã la c) pentru a = , b = , P(x) = şi Q(x) = x Deoarece z = + i este soluţie a ecuaţiei caracteristice λ + = , soluţia particularã se va cãuta sub forma x ( )cos ( )sin Py = xe ⎡⎣ Ax + B x + Cx + D x⎤⎦ Înlocuin în ecuaţie şi identificând coeficienţii se obţine sistemul A + D = , C = , − B + C = , − A = cu soluţia A = / , B = , C = , D = / Soluţia perticularã este deci −x cos x / + xsin x / Soluţia generalã a ecuaţiei este ( ) y x = C cos x + C sin x − x cos x / + xsin x / Ecua􀄠ii cu coeficien􀄠i variabili Pentru determinarea soluţiei generale a ecuaţiei omogene cu coeficienţi variabili nu exista metode generale Dacã însã aceastã soluţie poate fi precizatã, pentru determinarea unei soluţii particlare se poate folosi metoda variaţiei constantelor Teorema : Fie n n C y + C y + + C y solu􀄠ia generalã a ecua􀄠iei omogene ( ) ( ) ( ) ' n n n n y a x y a y ay − + + + − + = Dacã ( ) ( ) ( ) , , , n C x C x C x satisfac sistemul ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' n n n n n n n n n n n n n n C x y C x y C x y C x y C x y C x y C x y C x y C x y C x y C x y C x y f − − − − − − + + + = + + + = + + + = + + + = (x) ⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪ ⎩ atunci ( ) ( ) ( ) ( ) n n y x = C x y + C x y + + C x y este solu􀄠ie a ecua􀄠iei ( ) Exemplu : Sã se determine solu􀄠ia generalâ a ecua􀄠iei xy ''+ y ' = x, x > Se noteazã y ' = z Ecuaţia omogenã asociatã este xz '+ z = Rezultã z C x = adicã ( ) y = C ln x + C Se foloseşte metoda variaţiei constantelor pentru ( ) y x = ln x şi ( ) y x = Sistemul devine ( ) ( ) ' ln ' ' ' C x x C C x C x x ⎧ + = ⎪⎨ + ⋅ = ⎪⎩ cu soluţiile ( ) ( ) ln C x x C x x x ⎧ = ⎪⎨ = − ⎪⎩ Rezultã ( ) ( ) ln x C x A x x C x x B ⎧ = + ⎪⎪⎨⎪ = − + + ⎪⎩ Soluţia generalã a ecuaţiei este ( ) ln x y x = + A x + B Observa􀄠ie: Metoda varia􀄠iei constantelor poate fi folositã si pentru determinarea solu􀄠iilor particulare ale ecua􀄠iilor omogene cu coeficien􀄠i constan􀄠i atunci când func􀄠ia f (x) nu se încadreazã în situa􀄠iile prezentate anterior Exemplu: Sã se determine solu􀄠ia generalã a ecua􀄠iei '' ' ex y y y x − + = , x > Soluţia generalã a ecuaţiei omogene este ( ) x x G y x = C e + C xe Se aplicã variaţia constantelor pentru ( ) y x = ex şi ( ) y x = xex Sistemul obţinut este ( ) ' ' ' ' x x x x x x C e C xe e C e C e xe x ⎧ + = ⎪⎨ + + = ⎪⎩ cu soluţiile ( ) ( ) ln C x x A C x x B ⎧ = − + ⎪⎨ ⎩⎪ = + Soluţia generalã a ecuaţiei este y (x) = (−x + A)ex + xex (ln x + B) Ecua􀄠ii incomplete Ecuaţiile diferenţiale incomplete de ordin n sunt ecuaţii (în general neliniare) în expresia cãrora nu apar toate derivatele funcţiei necunoscute pânã la ordinul n- Penrtu rezolvarea lor se folosesc substitutii care le micşoreaza ordinul Ecua􀄠ii in care func􀄠ia necunoscutã apare doar prin derivatele sale Ecuaţiile de forma F(x, y(k ), y(k + ), , y(n) )= se reduc la o ecuaţie de ordin n-k prin substituţia z = y(k ) Exemplu : ( ) y''= + y' se transforma într-o ecuaţie de ordinul I folosind substituţia z = y' Ecuaţia z'= + z este o ecuaţie cu variabile separabile şi are soluţia generalã x ex C e C z = − − Rezultã deci cã ( ) ( ) x ex C e C y x z x dx C = = − − − ∫ Ecua􀄠ii în care variabila independentã nu apare explicit Aceste ecuaţii au forma F(y, y', , y(n) )= Se noteazã y'= p( y) Atunci ( ) ( ( )) ( ) ( ) ( ) ( ( ) ) ''' '' ' '' ' '' ' ' ' ' ' ' y y p p p p y y p y p y y p p y = = ⋅ + ⋅ = = = = ⋅ Si ecuaţia devine f (y, p, p', p(n− ) )= Aceeaşi procedurã se poate aplica pentru a micşora ordinul în continuare Sunt numeroase cazurile în care derivata apare în ecuaţie doar la puteri pare În acest caz se face notaţia ( ) p = y' Exemplu : ( ) y ⋅ y' '= + y' Se noteazã ( ) p = y' Rezultã cã p'(x) = p'(y(x))y'(x) = y'(x)y''(x) deci y''= p' / Din ecuaţia cu variabile separabile ( ) y p p '= + ⋅ rezultã p = Cy − adicã y'= ± Cy − Prin integrare directǎ de obţine C y Cy x A C y + ± = ⎟⎠ ⎞ ⎜⎝ ⎛ ln + − Soluţia ecuaţiei este ( ) ( Cx ) Cx C C e C e y ± ± − = Exerci􀄠ii propuse : Determina􀄠i solu􀄠iile generale ale urmãtoarelor ecua􀄠ii : y ''− y '+ y = x R : ( ) ( ) ( ) y x = C +C x e x + x + x + y ''− y '+ y = x + R : ( ) sin cos + + + ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ = + x x x C x y x e C x y ''+ y '+ y = e x R : ( ) ( ) y x = C + C x e−x + e x y ''− y '+ y = R : ( ) y x = C ex + C e x + y ''− y ' = ex R : ( ) x x x x y x = C e + C e− + e y ''+ y '− y = xe x R : ( ) x x x x y x = C e + C e− + x⎛⎜ − ⎞⎟e ⎝ ⎠ y ''+ y = cos x R : ( ) cos sin sin y x = C x + C x + x x y ''+ y = sin x R : ( ) ( ) cos sin cos y x = C x + C x + x y '''− y ''+ y = R : ( ) y x = C + C ex + C e x y '''− y = R : ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ + + = + − sin cos x C x e C y C e x x yIV + y = R : ( ) ( x x ) x x x C e C e y x C e C e − − + + = + + sin cos yIV − y '''+ y '' = ex R : ( ) x x y x C C x C C x e ⎛ ⎞ = + + ⎜ + + ⎟ ⎝ ⎠ '' ' e x y y y x − + + = , x > R : ( ) ( ) y x = C + xC e−x + xe−x ln x x y ''− xy '+ y = + x sin x R : ec omogenã are soluţia ( ) y x = C x + C x Ec iniţială are soluţia x x x x y C x C x cos sin − = + + − − x y ''+ xy '+ y = ln ( + x) R: ec omogenã are soluţia ( ) A B y x x x = + ( ) ( ) ln − − + − + = + + x x x x x C x C y CAPITOLUL Sisteme de ecua􀄠ii diferen􀄠iale de ordinul I Considera􀄠ii generale Forma explicitã a unui sistem de ecuaţii diferenţiale de ordinul I este ( ) ( ) ( ) ' , , , , ' , , , , ' , , , , n n n n n y f xy y y y f xy y y y f xy y y ⎧ = ⎪ = ⎪⎨⎪⎪ ⎩ = ( ) unde , , , n f f f sunt funcţii date, continue pe un domeniu din Rn+ O problemã Cauchy este formatã dintr-un sistem de ecuaţii diferenţiale şi un set de condiţii iniţiale, ( ) ( ) ( ) , , , n n y x = a y x = a y x = a ( ) O soluţie a sistemului ( ) este formatã din funcţiile , , , n y y y care verificã sistemul Astfel de sisteme de ecuaţii diferenţiale apar în mod natural în dinamicǎ, acolo unde legea fundamentalǎ ( F = m⋅ a ) se poate scrie (ţinând cont cǎ a = (x''(t ), y''(t ), z''(t ))şi ( ( , , , , ', ', '), ( , , , , ', ', '), ( , , , , ', ', ')) F = f t x y z x y z f t x y z x y z f t x y z x y z Ecuaţile de mişcare sunt deci ( ) ( ) ( ) ⎪⎩ ⎪⎨ ⎧ = = = '' , , , , ', ', ' '' , , , , ', ', ' '' , , , , ', ', ' z f t x y z x y z y f t x y z x y z x f t x y z x y z iar prin notaţia x'= u, y'= v, z'= w se obţine un system de ecuaţii diferenţiale de ordinul I In legaturǎ cu soluţile sistemelor de ecuaţii diferenţiale existǎ câteva rezultate importante Teorema Orice sistem de n ecuaţii diferenţiale diferenţiale de ordinul I scris sub forma ( ) este echivalent cu o ecuaţie diferenţialã de ordin n şi reciproc, orice ecuaţie diferenţialã de ordin n este echivalentã cu un sistem de ecuaţii diferenţiale de ordinul I Observa􀄠ie: soluţia sistemului ( ) poate fi gãsitã rezolvând ecuaţia de ordin n ataşatã prin metoda substituţiei ( se deriveazã una din ecuaţiile sistemului de n − ori, celelalte de n − ori şi se eliminã n − funcţii necunoscute) Aceastã metodã se mai numeşte şi metoda ecuaţiei rezolvante Exemplu : Sã se rezolve sistemul ' ' y z z y = ⎧⎨ ⎩ = − Derivând prima ecuaţie obţinem y '' = z ' Inlocuind aici z ' = −y (din a doua ecuaţie a sistemului) onţinem ecuaţia de ordinul II y ''+ y = cu soluţia ( ) y x = C cos x + C sin x Rezultã z(x) y'(x) C sin x C cos x = = − + Sisteme liniare şi omogene de ecua􀄠ii diferen􀄠iale cu coeficien􀄠i constan􀄠i Forma generalã a sistemului este ' ' ' n n n n n n n nn n y a y a y a y y a y a y a y y a y a y a y = + + + ⎧⎪ = + + + ⎪⎨⎪⎪ ⎩ = + + + ( ) Sistemului ( ) i se asociazã matricea coeficienţilor, anume ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎝ ⎛ = n nn n a a a a A Dacã notãm ( ) , , , n Y y y y τ = şi ( ) ' ', ', , ' n Y y y y τ = atunci sistemul se scrie în forma matricealã Y ' = A⋅Y ( ’) şi rezultatele prezentate la ecuaţii diferenţiale liniare de ordin n aratǎ cǎ mulţimea soluţiilor sistemului reprezintǎ un spaţiu vestorial de dimensiune « n » Soluţiile fundamentale ale sistemului vor fi cãutate sub forma ( ) i x i i i ni Y e= α α α τ λ unde i λ sunt valori proprii a matricii A, adicã soluţiile ecuaţiei caracteristice a sistemului : n n n n nn a a a a a a a a a λ λ λ − − = − ( ) iar constantele ij α trebuiesc determinate din sistem Dacǎ n Y , Y , ,Y sunt n soluţii liniar independente atunci soluţia generalã a sistemului este n n Y = C Y + C Y + + C Y ( ) Existã urmãtoarele situaţii importante : - Dacã ecuaţia caracteristicã ( ) are soluţiile reale şi distincte , , , n λ λ λ şi , , , n V V V sunt vectorii corespunzãtori acestor valori atunci soluţia sistemului este ( ) x x n x n n Y x CVe CV e CV eλ = λ + λ + + - Dacã ecuaţia caracteristicã ( ) are soluţii multiple (reale sau complexe), fiecare soluţie λ cu ordinul de multiplicitate p contribuie în suma ( ) cu termenii x Y Veλ = , ! x x Y = eλ ⎛⎜ V +V ⎞⎟ ⎝ ⎠ , …, ( ) ( )! ! ! p p x p p p x x x Y e V V V V p p λ − − − ⎛ ⎞ = ⎜⎜ + + + + ⎟⎟ ⎝ − − ⎠ unde , , p V V V sunt vestorii proprii corespunzãtori valorii proprii λ Problema rezolvãrii sistemului ( ) se reduce deci la determinarea valorilor proprii pentru matricea A şi a vectorilor proprii corespunzãtori acestor valori Metoda de rezolvare : a) Se scrie matricea A a sistemului ; b) Se determinã valorile proprii ale matricii rezolvand ecuaţia ( ) ; c) Pentru fiecare valoare proprie i λ se determinã vectorii proprii (atâţia cât e ordinul de multiplicitate al lui i λ şi se scriu soluţiile corespunzãtoare lui i λ ; d) Se scrie sistemul fundamental de soluţii al sistemului ; e) Se scrie soluţia generalã Exemple : Sã se determine solu􀄠ia generalã a sistemelor urmãtoare ' ' ' y y y y y y y y y y y y = − + ⎧⎪ = − + − ⎨⎪ ⎩ = − + a) Matricea sistemului este A ⎛ − ⎞ = ⎜ − − ⎟ ⎜ ⎟ ⎜ − ⎟ ⎝ ⎠ b) Valorile proprii ale matricii A sunt λ = , λ = şi λ = , adicã soluţiile ecuaţiei λ λ λ − − − − = − − c) Valoarea λ are ordinul de multiplicitate , deci va avea un singur vector propriu V = (α ,α ,α ) care verificã ecuatia α α α α α α ⎛ − ⎞ ⎛ ⎞ ⎛ ⎞ ⎜ − − ⎟ ⋅ ⎜ ⎟ = ⋅ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ − ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ ⎝ ⎠ ⎝ ⎠ Din rezolvarea sistemului compatibil nederminat cu un grad de libertate α α α α α α α α α α α α − + = ⎧⎪ − + − = ⎨⎪ ⎩ − + = ee obţine soluţia α α ⎛ ⎞ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ −⎝⎠ Se dã lui α o valoare particularã, de exemplu α = şi obţinem V ⎛ ⎞ = ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ −⎝⎠ In mod asemãnãtor obţinem V ⎛ ⎞ = ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ şi V ⎛ ⎞ = ⎜− ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ corespunzând valorilor λ şi λ d) Sistemul fundamental este x x x e Y e e ⎛ ⎞ ⎛ ⎞ ⎜ ⎟ ⎜ ⎟ = =⎜⎟ ⎜ ⎟ ⎜ − ⎟ ⎜ − ⎟ ⎝ ⎠ ⎝ ⎠ , x x x e Y e e ⎛ ⎞ ⎜ ⎟ =⎜ ⎟ ⎜⎜ ⎟⎟ ⎝ ⎠ , x x x e Y e e ⎛ ⎞ ⎜ ⎟ = ⎜ − ⎟ ⎜⎜ ⎟⎟ ⎝ ⎠ e) Soluţia generalã este C Y C Y C Y y y y = + + ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎝ ⎛ , adicǎ x x x x x x x x y Ce Ce Ce y Ce C e y Ce Ce Ce ⎧ = + + ⎪⎪ = − ⎨⎪ ⎩⎪ = + + ' ' ' y y y y y y y y y = + ⎧⎪ = + ⎨⎪ ⎩ = + Matricea sistemului este A ⎛ ⎞ = ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ Ecuaţia caracteristicã , λ − λ − = are rãdãcinile λ = şi λ =λ = − Un vector propriu al lui λ este V ⎛ ⎞ = ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ Subspaţiul valorii proprii λ =λ are dimensiunea Vectorii proprii satisfac ecuaţia α α α α α α ⎛ ⎞ ⎛ ⎞ ⎛ ⎞⎜ ⎟ ⋅⎜⎟ = − ⋅ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ ⎝ ⎠ ⎝ ⎠ , adicã sistemul cu douã grade de nedeterminare α α α α α α α α α + + = ⎧⎪ + + = ⎨⎪ ⎩ + + = Alegând α = − ,α = se obţine α = şi pentru α * = ,α * = se obţine α * = − Cei doi vectori proprii principali vor fi V ⎛− ⎞ = ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ⎠ şi V ⎛ ⎞ = ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ −⎝⎠ Valorile , , *, * α α α α se aleg astfel încât * * * * =α α −α α ≠ α α α α Soluţia generalã a sistemului este Y = CV e x +C V e−x +C (xV +V )e−x adicã ( ) x x x x x x x x y Ce Ce Cxe y Ce Ce y Ce Ce C x e − − − − − ⎧ = − − ⎪⎪ = + ⎨⎪ ⎩⎪ = + + − Sisteme liniare neomogene cu coeficien􀄠i constan􀄠i Forma generalã este Y '(x) = A⋅Y (x) + F (x) ( ) Ca şi în cazul ecuaţiilor liniare neomogene, soluţia generalã a sistemului neomogen este suma dintre soluţia generalã a sistemului omogen si o soluţie particularã a sistemului neomogen Pentru determinarea soluţiei particulare se poate folosi metoda variaţiei constantelor Teoremã : Dacã n n Y = C Y + +C Y este solu􀄠ia sistemului omogen asociat lui ( ) atunci o solu􀄠ie particularã a acestuia este ( ) ( ) P n n Y = C x Y + + C x Y unde func􀄠iile ( ) ( ) , , n C x C x satisfac ecua􀄠ia ( ) ( )( ) ( ) ( ) ' ' ' n n C x Y + C x Y + + C x Y = F x ( ) Din ecuaţia ( ) se calculeazã ( ) ( ) ' , , ' n C x C x şi apoi, prin integrare se obţin , , , n C C C Exemplu : S􀓽 se rezolve sistemul ⎪⎩ ⎪⎨ ⎧ = + + = − + −t t y x y e x x y e ' ' Sistemul omogen ⎩ ⎨ ⎧ = + = − y x y x x y ' ' are matricea coeficienţilor data de ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ − = A Valorile sale proprii sunt λ = şi λ = , vectorii proprii corespunzator acestora sunt ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ = V , respective ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ = V iar soluţia generalǎ a sistemului este ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ + + = ⋅ ⋅ + ⋅ ⋅ = ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ t t t t t t G G C e C e C e C e C V e C V e y x Expresiile pentru C şi C se determinǎ din sistemul ⎪⎩ ⎪⎨ ⎧ + = + = t t −t t t t C e C e e C e C e e ' ' ' ' Soluţii ale acestui sistem sunt / , / C = e− t − e−t C = −et − e− t O soluţie particularǎ a sistemului neomogen este ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − − − − = ⎟ ⎟⎠ ⎞ ⎜ ⎜⎝ ⎛ − − t t t t P P e e e e y x Rezultǎ cǎ soluţia generalǎ a sistemului este ⎪⎩ ⎪⎨ ⎧ = + − − = + − − − − t t t t t t t t y C e C e e e x C e C e e e Exerci􀄠ii propuse I Sã se rezolve sistemele urmãtoare : ' ' / y y z x z y z x + + = + ⎧⎪⎨ ⎩⎪ + − = R : ( ) y x = C e x +C e− x + x + x , z = −C e x + C e− x / − x / ' sin ' cos y y z x z y z x + + = ⎧⎨ ⎩ − − = R : sin ' sin cos z x y y y C ex C e x x x = − − = + − − − ' ' ' y z z u u y = ⎧⎪ = ⎨⎪ ⎩ = R : ( ) / / cos sin x x x x x y x = C e + C e− + C e− ( ) ( ) ( ) ( ) ' ' z x y x u x z x = = ' ' x x y y x y = − + ⎧⎨ ⎩ = − R : ( ) ( ) t t t t x t Ce Ce y t Ce C e − − = + = − ' ' x x y y x y t = + ⎧⎨ ⎩ = − + R : ( ) ( ) t t t t x t C C e t t y t C C e = + − − − = − + + − − ' ' ' y z u z y u u y z = + ⎧⎪ = + ⎨⎪ ⎩ = + R : ( ) ( ) ( ) ( ) ( ) ⎪ ⎪⎩ ⎪⎪⎨ ⎧ = + + − = + = − + − − − x x x x x x u x C e xC C C e z x C e C e y x C e C C e ⎪⎩ ⎪⎨ ⎧ = + = + − = − + z y z y x y z x x y z ' ' ' R : [ ( )] [ ( )] ( ) ⎪ ⎪⎩ ⎪⎪⎨ ⎧ = + + = + + + = + + t t t t t z C e C C t e y C e C C t e x C C t e ⎪ ⎪⎩ ⎪⎪⎨ ⎧ = − + − = − + − = − + / / / y y y y y y y y y y y y R : C e x C ex C e x y y y − ⋅ − ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎝ ⎛ − ⋅ + − ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎝ ⎛ ⋅ + ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎝ ⎛ − = ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎝ ⎛ ⎩ ⎨ ⎧ = − + = + − y x tgt x y tg t ' ' R : ⎩ ⎨ ⎧ = − + + = + + sin cos cos sin y C t C t x C t C t tgt ⎩ ⎨ ⎧ = + = − y x et x x y ' ' R : ( ) ( ) ⎪⎩ ⎪⎨ ⎧ = + + − = + + − t t y C C t t t e x C C t C t e II S􀓽 se rezolve problema Cauchy ⎪ ⎪ ⎪ ⎩ ⎪ ⎪ ⎪ ⎨ ⎧ = = − + ⋅ = ⋅ = + − ( ) , ( ) / / y z z x x z y y x x y R : ( ) ( ) ⎪⎩ ⎪⎨ ⎧ = + − = + x x z x x y x x CAPITOLUL Elemente de calcul opera􀄠ional şi aplica􀄠ii în teoria ecua􀄠iilor diferen􀄠iale Calculul operaţional se ocupã cu studiul transformãrilor integrale Acestea, numite şi “operatori integrali”, transformã derivarea şi integrarea în operaţii algebrice Ecuaţiilor diferenţiale şi integrale le corespund ecuaţii algebrice Pentru a rezolva o ecuaţie diferenţialã este suficient sã se rezolve ecuaţia algebricã şi sã se aplice transformarea integrală inversă soluţiei obţinute Cele mai directe aplicaţii în studiul ecuaţiilor diferenţiale îl are transformata Laplace Transformata Laplace Transformata Laplace este un operator între douã spaţii de funcţii, operator care transformã derivarea şi integrarea în operaţii algebrice Notãm FR ={ f : R→R} Defini􀄠ie: Func􀄠ia R f ∈F este un original dacã satisface condi􀄠iile urmãtoare: a) f (t ) = pentru orice t şi s > astfel încât | f (t ) | Numãrul pozitiv { ( ) } s = min s | | f t | se numeşte indice de creştere (sau abscisã de convergen􀄠ã) Mul􀄠imea func􀄠iilor original se noteazã R O Exemple: Urmãtoarele funcţii sunt funcţii original: a) ( ) , , ekt t f t t ⎧⎪ ≥ = ⎨ ⎩⎪ ⎧⎪ = = ⎨⎪ ⎩ Ω soluţia generalǎ este ( ) t t g x t C e C e ⎟⎠ ⎞ ⎜⎝ ⎛ Ω − − − ⎟⎠ ⎞ ⎜⎝ ⎛ − + −Ω = + α α α α - dacǎ α = Ω soluţia generalǎ este ( ) ( ) t g x t = C t +C e−α - dacǎ α aratǎ cǎ lim ( ) = →∞ x t g t adicǎ, în absenţa unui termen perturbator, mişcarea este amortizatǎ Deoarece ecuaţia caracteristicǎ nu are rǎdǎcini pur imaginare, soluţia particularǎ se cautǎ sub forma x (t ) A t B t p = sinω + cosω Inlocuind x'(t ) = Aω cosωt − Bω sinωt şi x''(t ) = −Aω sinωt − Bω cosωt în ecuaţia iniţialǎ obţinem, prin identificarea coeficienţilor, ecuaţiile ( ) ⎪⎩ ⎪⎨ ⎧ + Ω − = Ω − − = ( ) A B A B a αω ω ω αω din care rezultǎ ( ) ω α ω ω Ω − + = Ω − A a şi ( ) ω α ω αω Ω − + B = −a In acest caz soluţia particularǎ este ( ) ( ) ( ) x t a t a t p ω ω α ω ω αω ω α ω ω sin cos Ω − + − Ω − + = Ω − Ca de obicei, soluţia generalǎ a ecuaţiei este x(t ) x (t ) x (t) g p = + Exemplu: Sǎ se resolve ecuaţia x"+ x'+ x = sin t In acest caz α = , Ω = , ω = Ecuaţia caracteristicǎ r + r + = are soluţiile , r = − ± i , deci soluţia generalǎ a ecuaţiei omogene este x (t ) e t (C t C t ) g cos sin = − + Soluţia particularǎ este x (t ) t t p cos sin = − Ea se mai poate scrie sub forma x (t ) = ( t −φ ) p sin , unde tgφ = Ea are perioada principalǎ π π ω = π = = T care este consideratǎ perioada de oscilaţie Sǎ observǎm însǎ cǎ soluţia generalǎ x(t ) e t (C t C t ) ( sin t cos t ) cos sin = − + + − Nu este o funcţie periodicǎ Mişcarea unui pendul (ecua􀄠ie liniar􀓽) S􀓽 se descrie mişcarea unui pendul legat de punctul O printr-un fir de lungime l, dac􀓽 acest pendul este deviat de la verticala locului cu unghiul la centru θ şi viteza sa ini􀄠iala este v Figura Pendulul matematic Se considerǎ ca variabilǎ independentǎ unghiul dintre verticala locului şi firul OM (pendulul se aflǎ în punctul M ) Funcţia necunoscutǎ x = x(θ ) reprezintǎ arcul OM Asupra corpului aflat în M acţioneazǎ greutatea G = mg dirijatǎ pe verticalǎ în jos şi tensiunea în fir dirijatǎ de-a lungul firului spre punctul O Forţa rezultantǎ este F = −mg sinθ ⋅ t + (T − mg cosθ )⋅ n , unde t, n reprezintǎ vectorul tangent, respectiv normal, la traiectorie Proiectând legea fundamentalǎ F = m⋅ a de-a lungul tangentei la traiectorie şi tinând cont cǎ acceleraţia este a doua derivatǎ a spaţiului, obţinem ecuaţia θ θ M T A O P x''(θ ) = −g sinθ Dacǎ θ (adicǎα > Ω ) atunci soluţia generalǎ este ( ) t t q t A e Be ⎟⎠ ⎞ ⎜⎝ ⎛ − Ω − − ⎟⎠ ⎞ ⎜⎝ ⎛ −Ω − = ⋅ + α α α α Constantele A şi B se determinǎ din condiţiile iniţiale ( ) q = q şi i( ) = Din sistemul ⎪⎩ ⎪⎨ ⎧ = ⎟⎠ ⎞ ⎜⎝ ⎛ ⎟⎠ ⎞ ⎜⎝ ⎛ + Ω − − ⎟⎠ ⎞ ⎜⎝ ⎛ −Ω − + = A α α B α α A B q rezultǎ −Ω = ⋅ −Ω + α q α α A şi −Ω = −Ω − α q α α B Sǎ remarcǎm cǎ puterile ce apar la exponenţi sunt amândouǎ negative şi deci lim ( ) = →∞ q t t - dacǎ C L R = ⋅ (adicǎα = Ω ) atunci soluţia generalǎ este q(t ) = (At + B)⋅ e−αt Coeficienţii A şi B sunt în acest caz A = −q şi B = q şi soluţia se scrie q(t ) = q (t − )⋅ e−αt - dacǎ C L R Ecuaţia curbei cǎutate va fi y = y(x) Raza de curburǎ este definitǎ de ( ( ) ) '' ' / y ρ = − y , deci ecuaţia problemei este ( ( ') ) '' / − y = R⋅ y Aceastǎ ecuaţie neliniarǎ incompletǎ se transformǎ într-un sistem de ecuaţii de ordinul I ( ) ⎪⎩ ⎪⎨ ⎧ ⋅ = + = / ' ' R u u y u Din ultima ecuaţie a sistemului rezultǎ, prin integrare directǎ R x C u u = − + adicǎ R (x C) x C u − − = ± − Inlocuind u în prima ecuaţie şi integrand obţinem ( ) y = − R − x −C +C , adicǎ ( ) ( ) y −C + x −C = R Aceastǎ relaţie aratǎ ca doar cercurile de razǎ R (şi centru arbitrar ( ) A C,C ) au raza de curburǎ constantǎ, egalǎ cu R Bibliografie [ ] Bălan T , Transformata Laplace, Ed Universitaria, Craiova, [ ] Brînzănescu V , Stănăşilă O , Matematici Speciale – teorie, exemple, aplicaţii, Ed All, Bucureşti, [ ] Constantinescu D , Équations Différentielles, Chapitres de Mathématiques Supérieures, Ed Universitaria, [ ] Corduneanu A , Ecuaţii diferenţiale cu aplicaţii în electrotehnică, Ed Facle, Timişoara, [ ] Craiu M , Roşculeţ M , Ecuaţii Diferenţiale Aplicative, EDP, Bucureşti, [ ] Matei Stefan, Analiza circuitelor, Editura Didactică şi Pedagogică, [ ] Micula Gh , Paraschiva P , Ecuaţii diferenţiale şi integrale prin probleme şi exerciţii, Ed Dacia, Cluj-Napoca, [ ] Olariu V , Stănăşilă O , Ecuaţii diferenţiale şi cu derivate parţiale, Ed Tehnică, Bucureşti, [ ] Predoi M , Constantinescu D , Racilă M , Teme de calcul diferenţial / Teme de calcul integral, Ed Sitech, Craiova, [ ] Rădoi M , Deciu E , Mecanică, Ed Didactică şi Pedagogică, Bucureşti, [ ] Roşca D , Mecanică, Ed Didactică şi Pedagogică, Bucureşti, [ ] Săvescu M , Metode de analiză a circuitelor electronice, Ed Stiinţifică şi Enciclopedică, Bucureşti, [ ] Soare M , Teodorescu P P , Toma I , Ecuaţii diferenţialw cu aplicaţii în mecanica construcţiilor, Ed Tehnică, [ ] Turcitu G , Şterbeţi C , Matematici Speciale – Analiză complexă şi ecuaţii diferenţiale, Ed Radical, Craiova, 